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abstract 

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts 
KJi , and observed with the Herschel Space Observatory as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered 

■ Extragalactic Survey (HerMES) guaranteed time key programmes. Our study makes use of the broad far-infrared and submillimetre 

■ wavelength coverage ( 100-600 yum) only made possible by the combination of observations from the PACS and SPIRE instruments 
aboard the Herschel Space Observatory. Using a power-law temperature distribution model to derive infrared luminosities and dust 
temperatures, we measure a dust emissivity spectral index for SMGs ofp = 2.0 ± 0.2. Our results unambiguously unveil the diversity 
of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~ IO'^Lq and relatively warm dust components, while 
others are fainter (a few times IO'-Lq) and are biased towards cold dust temperatures. Although at z ~ 2 classical SMGs (> 5 mjy 
at 850/im) have large infrared luminosities (~ lO'^L©), objects only selected on their submm flux densities (without any redshift 
informations) probe a large range in dust temperatures and infrared luminosities. The extreme infrared luminosities of some SMGs 
(Lm ^ 10'--^ Lq, 26/61 systems) imply star formation rates (SFRs) of > 500 M© yr"' (assuming a Chabrier IMF and no dominant AGN 
contribution to the FIR luminosity). Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead 
correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is 
the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with 
observations of local ultra-luminous infrared galaxies triggered by major mergers and with results from hydrodynamic simulations 
of major mergers combined with radiative transfer calculations. Moreover, we find that luminous SMGs are systematically offset 
from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty 
cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have 
cold dust temperatures, are located close to the main sequence of star formation, and therefore might be evolving through a secular 
mode of star formation. However, the properties of this latter population, especially their dust temperature, should be treated with 
caution because at these luminosities SMGs are not a representative sample of the entire star-forming galaxy population. 
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1. Introduction 

Submillimetre (submm) observations probe the Rayleigh-Jeans 
side of the blackbody emission of dust in galaxies. In that 
regime, the dimming of the submm flux density of a galaxy 
due to its cosmological distance is counterbalanced by the red- 
shifting of its spectral energy distribution (SED). Consequently, 
submm observations can trace galaxies with the same infrared 
luminosities over a broad range of redshifts, and are thus a very 
powerful too l for studying the cosmic star-formation history 
(lBlainllT996l) . Unfortunately, most current deep submm surveys 
have spatial resolutions on the order of ten arcseconds. This 
large b eam size, combined with the steep submm number counts 
(e.g., ICoppin et al. 2006^, leads to a high level of confusion, 
which ultimately limits the sensitivity of submm observations. 
Submm surveys are therefore limited to th e brightest s ources 
and sub mm-selected galaxie£] (SMGs; 'Smail et alJ Il997t 
iBarger et al. 1998; Hughes et al. 1998; Blain et al. 200^ for a 
review) have thus been primarily used for probing the most lu- 
minous tail of the high-redshift star-forming galaxy population. 

Substantial efforts have been invested in high-resolution 
multi-wavelength identifications of SMGs using (sub)mm, 
radio, mid- or near-infrared ob servations (e.g., Downes et al. 
1999'; 'Dannerbaueretal.' '2002'; Ivison et al. 2002t 

2005: Bertoldi et al.. 2007: Biggs et al.. .20lC ~It 

found that SMGs lie at high-redshift, z ~ 2 (iHughesetal 



IPope et alj 
has been 



1998 



2005; 



'Barger et al.' '2000'; 'Smail et al 
Pope et al. 2006; Wardlow et_al 



Carilli & Yun 
2000[ IChapman et al 

201 ih . and ar e massive systems (M, ~ 10^"-10'^ 



Swinbank et al]|2004l; iTacconi et alj|2006l 120081; iHainUne et all 
201 ih . Extrapolation of their infrared luminosities (Lir) from 
submm, radio or mid-infrared observations, have shown that 
SMGs are extremely lumin ous (Lir(8 - 1000 jum) > 10'^ Lq; 
e.g.. Chapman et al. 2005; Pope et alJl2006l l2008t iKovacs et alJ 
l2006i ,20 10) . Their infrared luminosities are mainly powered 
by star-f ormation rather than by active galactic nucleus (AGN) 
activi ty ('Alexander et al.' 2005; Lutz et al. 2005; Valiante et alj 
20071; rMenendez-Delmestre et al. 2007; Pope et al. 200l 
Menendez-Delmestre et al. 2009; Laird et al. 2010), and corre- 
spond to star-formation rates (SFRs) of a few 100s to few 1000s 
of Moyr"'. The most luminous SMGs are therefore peculiar 
galaxies because their SFRs are higher than that of typica l 
galaxies of similar mass at similar redshift (iDaddi et alJl2007bl) . 
Interferometric observations of their CO molecular gas suggest 
that the most luminous z ~ 2 SMGs (flux density at 850 fim, 
5 850 > 5 mJy) are major mergers in various stages, characterised 
by compact or very distu r bed CO kin ematics/morph ologies 



by compact or very distu r bed CU Kin ematics/morph ologies 
(ITacconi et al.1 12006[ 120081; lEngel et alj [2010: Bothwell et alj 
120101) . The gas to total baryonic mass fraction of SMGs is 
comparable to that of typical galaxies at the same redshift 
(30 - 60%; Tacconi et al. 2008, 2010), impjying that SMGs 
have higher star-fo rmation efficiencies dPaddi et al.ll2008ll201(i: 
iGenzel et alj 120101) . Finally, although the comoving volume 
density of SMGs wi th Ssso > 5 mJy is low (~ lO^^Mpc"^; 
IChapman et al.ll2005l) . their contribution to the SFR density of 
the Universe at z ~ 2 is ~ 10% (Chapman et al. 2005). 

Based on these derived properties, a picture of the na- 



Send offprint requests to: B. Magnelli, e-mail: magnelli@mpe.mpg . de 
* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

' Note that here we use the term SMGs to refer to sources selected 
by ground-based facilities in the 850-1200 yum window. 



ture of the most luminous SMGs has emerged. SMGs with 
S 850 > 5 mJy are thought to exhibit very intense short-lived 
star-formation bursts, triggered by mergers, and to be the 
high-redshift prog enitors of local massi ve early-type galaxi es 
(Lillv et al.1 Il999li ISwinbank et alj " 
Tacconi et al. 2008; Cimatti et all 



IDaddi et all iBoTallbl; 
In fliat picture, SMGs 
belong to a class of galaxies offset from the so-called "main 
sequence of star-formation" which links the SFRs and stellar 
masses of normal st ar-forming ga l axies ( SFGs) over abroad 



ran ge of r edshifts ("Noeske et al.1 120071: iDaddi et all l2007bt 



Elb az et al.l l%07; Rodighiero et aLl l2010ri201 Ih . The existence 
of this main sequence of star-formation is usually interpreted 
as evidence that the bulk of the SFG population is forming 
stars gradually with a long duty cycle, likely sustained by the 
accretion of cold gas fro m the interga lactic medium (IGM) and 
along the cosmic web dPekel et alJ 2009; Dave et al. 2010|). 
Occasional major merger events create extreme systems with 
intense short-lived starbursts, like SMGs, which are offset from 
the main sequence of star-formation and which likely evolve 
into "red and dead" galaxies. 

The picture of SMGs as a homogeneous population of 
major mergers has now been weakened by new observational 
constraints. The (sub)mm selection method does not correspond 
to a perfect bolometric selection but rather selects galaxies in 
the Tdust - LiYi parameter space favouring, at lo w infrared lumi- 
nosities, galaxies with cold er dust temperature (IChapman et al.l 
120051: iMagnenietaLl 120101) . Thus, current SMG samples can 
contain a significant fraction of relatively low luminosity galax- 
ies with cold dust temperature, i.e., galaxies with lower SFRs in 
the main sequence regime. The diversity of the SMG population 
is also supported by high-resolution observations. Some submm 
sources are actually composed of two galaxies (with normal 
ongoing star-formation) which are soon to merge and are 
observ ed as o ne sub mm so u rce because of the l a rge submm 
beani dYounger et al.l 120091; |Kovacs et alj 120101; IWang et all 
1201 ll) . Finally, constraints from simulations also support this 
diversity. While simulations of major mer gers are able to 
reproduce the extreme SFRs of b right SMGs ([Chakrabarti et al.l 
2008; Naravar ianet al.| [2010; Ha vward et al.ll2011h . fliere might 
be issues (depending on the exact merger condition needed to 
create these properties) to match the comoving volu me density 
of SM Gs usin g the high-r e dshift major merger rates (iDave et al.l 
12010 1). Thus, iDave et al ] (120101) have tried to reproduce the 
properties of SMGs using hydrodynamic simulations in a 
cosmological context. Their simulations cannot simultaneously 
reproduce the measured SFRs and comoving densities of SMGs, 
because the bulk of their simulated SMGs evolve secularly and 
exhibit lower SFRs than those infeiTed from observations (by a 
factor ~ 2 - 3). These results are also consistent with those of 
semi-analytic models which have great difficulties accounting 
simultaneously f or the measured luminosities/SFRs and nu mber 
counts of SMGs dBaugh et alJl2005LISwinbank et al.ll2008l) . 

Due to all these difficulties some questions remain: How 
homogenous is the SMG population? Have SMG luminosities 
been overestimated? What triggers their SFRs? 

One of the ingredients needed to shed light on the nature 
of SMGs is direct and robust measurements of their infrared 
luminosities and SEDs. Indeed, while SMGs have been studied 
at all wavelengths, in most cases their infrared luminosities are 
still based on large extrapolations from radio, submm or mid- 
infr ared obse r vation s. Using 350 yum SHARC-2 observations, 
Kov acs et al. I d2006l UOIO) provided more robust estimates 
of the infrared luminosity of a handful of SMGs. However, 
these studies still lacked rest-frame far-infrared observations on 
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both sides of the peak of the SEDs. Using observations by the 
1.8-m Balloon-borne Large Aperture Submillimetre Telescope 
(BLAST) at 250, 350, 500 jum, Chapinetal. (2011) studied 
the far-infrared SED of SMGs at its peak and thus robustly 
constrained their dust temperatures. Nevertheless, this study 
was limited to a relatively small SMG sample (23 sources with 
spectroscopic redshift estimates) and suffered from observations 
with large beam size (i.e., ~ 19" at 250//m). N ow, thanks to 
the a dvent of the Herschel Space Observatory dPilbratt et alj 
we can go further in the analysis of the far-infrared 
SED of SMGs. Using deep observations at 100 and 160 //m 
by the Ph otodetector Array Camera and Spectrometer (PACS; 
Poglitsch et al. 2010) onboard the Herschel Space Observatory, 
Magnelli^eLaD (12010|) estimated the infrared luminosities and 
dust temperatures of a small sample of SMGs (17 sources). Soon 
after, [Chapman et al. (2010) provided similar estimates using 
deep observations at 250, 350 and 500 /vm using the Spectral 
and Photometric Imaging REceiver (SPIRE; Griffin et al...201ft) 
also on Herschel. Both studies revealed the diversity of the 
SMG population and its bias, with respect to a bolometric 
selection, towards galaxies with cold dust temperature. Some 
galaxies exhibit extreme infrared luminosities of ~ lO'-'L© 
and relatively warmer dust components, while others have 
much lower luminosities (i.e., a few 10'^ Lq) and colder dust 
components. 

After more than two years of operation, Herschel has now 
produced deep observations of the most widely studied blank 
and lensed extragalactic fields. These combined new PACS and 
SPIRE data provide for the first time a wide coverage of the 
far-infrared SEDs of a large sample of SMGs, allowing us to 
go further in our understanding of their properties. Our results 
unambiguously reveal the true infrared luminosity of SMGs and 
can be used to test the quality of pie-Herschel estimates based 
on monochromatic extrapolations. These infrared luminosities 
and dust temperatures also shed light on the diversity of this 
population and can be used to test the different modes of star 
formation that could power their luminosities. Finally using the 
large wavelength coverage provided by the Herschel observa- 
tions, we can constrain the dust emissivity spectral index, /3, of 
SMGs. 

Here, we use PACS and SPIRE data for a sample of 
61 SMGs with known spectroscopic redshifts to provide an 
insight into the properties and nature of the SMG population. 
A comprehensive analysis of the complete SMG samples in the 
fields studied here will be the subject of other papers. 

The paper is structured as follows. In section |2] we present 
the Herschel data used in our study. Section [3] presents our 
He rschel-detected SMG sample with known spectroscopic 
redshifts and discusses the selection function of this sample. 
Section |4] is dedicated to SED analysis, describing how we 
have derived dust temperatures and infrared luminosities using 
a single-temperature modified blackbody model and a power- 
law temperature distribution model. We consistently refer to 
temperatures as T^asi if based on a y6 - 1.5 modified black- 
body, and Tc for the minimum temperature in the power-law 
distribution model. Scientific conclusions drawn from these 
estimates are discussed in Section|5]and in section|6]we discuss 
the nature of SMGs. Finally, we summarize our findings in 
Section |7] Throughout the paper we use a cosmol ogy with 
Hp ^ 7Ikms"'Mpc"', Qa = 0.73 and Qm = 0.27. A lChabried 
(I2003h initial mass function (IMF) is always assumed. 



2. Observations 

In this study, we used deep PACS 70, 100 and 160jum and 
SPIRE 250, 350 and 500 jum observations provided by the 
Herschel Space Observatory. PACS observation s were taken as 
part of the PACS Evolutionary Probe (PElfl Lutz et alJl20TTI) 
guaranteed time key programme, while the SPIRE observations 
were taken as part of the Herschel Multi-tiered Extragalactic 
Survey (HerME^J; Oliver et al. 2011, MNRAS submitted.). 
These two large key programmes are structured as "wedding 
cakes" (i.e., with large area wide surveys and smaller pencil 
beam deep surveys) and include many widely studied blank and 
lensed extragalactic fields. Many of these fields being common 
to both programmes, their combination provides an unique and 
powerful tool to study the SED of galaxies over a broad range of 
wavelength. The PEP and HerMES s u rveys and data reduction 
methods are described in ' Lutz et aP (1201 ih and Oliver et al. 
(201 1, MNRAS submitted) and references therein, respectively. 
Here, we only summarise the properties relevant for our study. 

From the PEP and HerMES programmes, we used 
the observations of the Great Observatories Origins Deep 
Survey-North (GOODS-N) and -South (GOODS-S) fields, 
the Lockman Hole (LH) field, the Cosmological evolution 
survey (COSMOS) field and the lensed fields Abell 2218, 
Abell 1835, Abell 2219, Abell 2390, Abell 370, Abell 1689, 
MS 1054, CL0024 and MS045. Table [1] summarises the main 
properties of these fields. Herschel flux densities were derived 
with a point-spread-function-fitting analysis guided using the 
position of sources detected in deep 24 fim observations from 
the Multiband Imaging Photometer (MIPS: iRieke et al.l l2004l) 
onboard the Spitzer Space Observatory. This method has the 
advantage that it deals with a large part of the blending issues 
encountered in dense fields and providing a straightforward 
association between MIPS, PACS and SPIRE sources. This 
MIPS-24jum-guided extraction is also very reliable for the 
purpose of this study, because here we focus on a subsample 
of SMGs which al ready have, for the m ost part, a MIPS-24yL(m 
identification (e.g.. lHainUne et al]|2009l) . 

In PEP, prior source extraction was performed using the 
method presented in Magnelli et al. ( 2009), wh ile in HerMES it 
was performed using the method presented in iRoseboom et al.l 
(2010), both consortia using consistent MIPS-24 jum catalogues. 
In GOODS-N and -S, we used die GOODS MI PS-24 yum 
catalogue presented in iMagnelli et al.l (l2009l I20TI reaching 
a 3cr limit of 20yuJy. In the LH, we used the MIPS-24 ;um 
catalogue provided by a Spitzer legacy programme (PI: E. 
Egami), reaching a 3cr limit of 30 juJy (Egami et al., in prep.). In 
COSMOS, we used the latest M IPS-24 jUm catalogue a vailable, 
reaching a 3cr limit of 45 /iJy (iLe Floc'h et all |2009'). In the 
lensed fields, we used the public MIPS-24 /im observations 
(PI: G. Rieke). The data processing and catalogue extraction 
follow the standard MIPS processing with some improvements, 
this is described in more detail in Valtchanov et al. (in prep.). 
In the central region these MIPS-24 //m data reaches a Icr 
limit of ~ 20-100;uJy depending on t he cir rus contamination 
(e.g., iMarcillac et akl 120071: iBai et al.1 l2007h . Using all these 
MIPS-24 yum source positions as prior, we created our PACS 
and SPIRE catalogues. The reliability, completeness and con- 
tamination of our PACS and SPIRE cata logues were tested via 
Monte-Carlo simulations (see Lutz et al.. l20I ll and Oliver et al. 



^ http://www.mpe.mpg.de/ir/Research/PEP 
' http://hermes.sussex.ac.uk 
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2011, MNRAS submitted for details). All these prop erties are 
given in lBerta et af] (1201 ih and'Rose boom et alJ (l2010t) . Table [l] 
only summarises the depth of all these catalogs. 

We note that the SPIRE prior catalogues reach a 3cr limit 
of ~ lOmJy, ~ 12mJy and ~ 15mJy at 250, 350 and 500 yum, 
respectively, while the formal 3cr extragalactic confusion limits 
at these wav elengths are 14.4 mJy, 16.5 mJy and 18.3 mJy 
dNguven et al. 2010). Sources detected below these formal 
3cr confusion limits should thus be treated with caution. In 
our specific SMG sample (with robust spectroscopic redshift 
estimates), only a small fraction of galaxies has SPIRE mea- 
surements below these formal confusion limits (less tha n 10%). 
For th ese sources, we follow the prescription of Elbaz et alJ 
(l20Toh . i.e., we take advantage of the higher spatial resolution 
of the MIPS-24yum observations to flag some galaxies as more 
"isolated" than others and for which SPIRE flux densities can 
potentially be more robust. Using this diagnostic, we conclude 
that in our final SMG sample only three sources (i.e., 5% of 
our sample) have SPIRE measurements potentially affected by 
confusion. While useful, we note that this diagnostic might 
not be fully reliable in fields where only shallow MlPS-24yum 
observations are available. In our case, only the COSMOS field 
can significantly be affected by this limitation and in this field 
none of our SMGs only relies on SPIRE flux densities below 
the formal 3cr confusion limit. 



3. Galaxy sample 

In order to infer dust temperatures, infrared luminosities and 
more generally dust properties, we have to rely on SMGs 
with robust redshift estimates obtained through secure multi- 
wavelength identifications. In this section, we present the con- 
struction of such a sample and discuss its selection function. 

In every field the construction of our sample follows three 
steps, (i) First, we search in the literature for samples of SMGs, 
i.e., galaxies selected by ground-based facilities in the 850- 
1200;um window, with robust multi-wavelength identifications 
and spectroscopic redshift estimates. In some of our fields, 
more than one such SMG sample were available. For example 
in GOODS-N, multi-wavelength identificat ion of Submillimetre 
Common User Bolometer Array (SCUBA: [Holland et aljri999l) 
and AzTEC (Wilson et al. 2008) sources have been separately 
published. In that case, we cross-match these samples using a 
matching radius of 9" (i.e., about the half-width at half max- 
imum, HWHM, of the submm observation^ and keep, for 
sources presented in more than one sample, the more secure 
multi-wavelength identifications (i.e., the one with the lo west 
probability, P, of chance association. iDownes et al.iri986h . (ii) 
We complement the far-infrared SED coverage of the SMGs de- 
fined in step (i) by searching for their submm/mm counterparts 
in all blind catalogues available (i.e., catalogues with no multi- 
wavelength identifications). In this step we again use a match- 
ing radius of 9". (iii) Finally, we cross-match the SMG sam- 
ple defined in step (i) (and which SED coverage has been com- 
plemented in step (ii)) with our MIPS-PACS-SPIRE catalogues. 
In this step we use the optical, MIPS or radio positions of the 
SMGs, the MlPS-24 fim positions from our MlPS-PACS-SPlRE 
catalogues and a matching radius of 3" (i.e., corresponding to 
the MIPS-24 fim HWHM). 

Some of our SMGs with robust spectroscopic redshift esti- 

This radius also corresponds to the 3crpos positional error of submm 
observations (cTpos ~FWHM/(2 x SN)), assuming that the bulk of our 
submm detections has a signal to noise ratio (SN) of ~ 3. 



mates might correspond to a PACS/SPIRE detection missed by 
our source extraction method because of a lack of a MlPS-24 /im 
prior. For that reason, we visually check in our PACS/SPIRE im- 
ages that the absence of a PACS/SPIRE detection was not due to 
a lack of a MlPS-24 fim prior. We find no such cases. 

3.1. GOODS-N 

In GOODS-N, we use the multi-w avelength i d entification o f 
SCUBA-850//m sources made by ' Pope et all (|2006' '20080 
using data a nd redshift info r mation s mainly from^orvs et aP 
(120031) and IChapman et aP (l2005l) . We also use the multi- 
wavelength i dentification of AzTEC- 1.1 mm sources made by 
IChapin et al.l (2009). From the Pope et al. sample we only use 
the SMGs with spectroscopic redshift estimates. 

From the AzTEC sample of Chapin et al., we only keep 
the two sources with robust spectroscopic redshifts that are not 
detected by SCUBA (i.e., not already included in the Pope et al. 
sample). We complement the Pope et al. sample with AzTEC 
flux densities wh en available. 

iGreve et all ( [2008) present the Max Planck Milhmeter 
Bolometer (MAMBO, at 1.2 mm) observations of the GOODS- 
N field. Some of these MAMBO sources have robust radio 
identifications in this paper but the corresponding radio positions 
are not provided. Consequently we only consider MAMBO 
counterparts of our SCUBA and AzTEC sources. 

This sample of 25 SMGs with robust redshift estimates is 
cross-matched with our MIPS-PACS-SPIRE multi-wavelength 
catalogue. Fourteen SMGs are detected in at least one of the 
PACS-SPIRE bands. Among those 14 sources, 10 are detected 
by both PACS and SPIRE, 3 are only detected with SPIRE and 
1 only detected with PACS. The final sample of 14 SMGs in 
GOODS-N is presented in Tables[3]and|4l 

Four SMGs are detected only in the SPlRE-250 /vm band 
with flux density below the formal 3cr confusion limit, namely, 
GN5, GN15, GN20 and GN20.2. For these four sources we 
compute their "cleanness" index as defined in lElbaz et aP 
(2010), i.e., sources are defined as "isolated" if they have at 
most one MlPS-24 //m neighbour within 20" with ^24 > 50% 
of the central MlPS-24 yum source. Among those four sources, 
one is found to be "isolated" and hence with robust SPIRE 
measurements (GN15). Therefore, results derived for GN5, 
GN20 and GN20.2 have to be treated with caution. 

3.2. GOODS-S 

In GOODS-S we use the multi-wavelength identification of 
sources observed by the Large APEX Bolometer Camera 
(LABOCA) ECDFS Submm Survey at 870 //m (LESS; 
Wei B et al.1 12009b). as presented by Big gs et alJ dlOll). This 
sample contains 75 SMGs robustly associated to MIPS, radio 
and optical counterparts but only 15 are situated in the deep 
GOODS-S field observed by Hersche^ Redshift information is 
taken from zLESS (Danielson et al., in prep.) which provides 
spectroscopic follow-up of the WeiB et al. sources. 



' For GN05 we use the spectroscopic redshift revised in iPope et al.l 
(2008); for GN20 and GN20.2 we use the spectroscopic redshifts re- 
vised in Daddi et al. (2009a b); and finally for GN07 we use the redshift 
from lChapman et al. (2005) 

* PEP and HerMES have both observed the Extended Chandra Deep 
Field South. These observations are shallower than those of GOODS-S 
and are not used in this analysis. 
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IScott et alJ (1201 oh presented the AzTEC observations of 
the GOODS-S field, but no multi-wavelength identifications of 
these sources are available. 

This yielded seven SMGs with robust spectroscopic 
redshift estimates. This sample is then cross-matched with 
our MIPS-PACS-SPIRE multi-wavelength catalogue. These 
seven SMGs are all detected in at least one PACS/SPIRE band. 
Multi-wavelength properties of these seven SMGs are presented 
in Tables |5] and |6] 



3.3. Lockman Hole (LH) 

In LH, we start from the multi-wavelength identifications of 
44 SCUBA HAlf Degree Extragalactic Survey (SHADES ; 
ICoppin et al.l 12006 *) sources made by Ivison et al.' (2001). 
Eleven have a spectroscopic redshift in Chapman et al. (2005). 
These SCUBA sources were associated in Ivison et al. (120071) 
with their MAMBO counterparts (Greve et al. 200 4). We also 
used the AzTEC counterparts of these sources provided in 
lAust ermann et al.l (1201 Oh. 

Chapman et aij( 2005h provide redshift information for two 
additional SCUBA SMGs that are not in the Ivison et al. sample 
(SMMJ105225.79-H571906.4 and SMMJ105238. 19-1-571651.1). 
The absence of these two SMGs in this sample could be 
explained by their low S/N submm detections. We decided to 
include those two galaxies in our sample of SMGs with robust 
redshift estimates; 

Recently, ICo ppin et alj (1201 Ol) derived the spectroscopic 
redshifts of six SMGs using the PAH signatures observed 
in the Spitzer-IRS spectra. This study added one SHADES 
source (LOCK850.15) and four AzTEC sources (AzTEC.Ol, 
AzTEC.05, AzTEC. 10 and AzTEC.62) to our SMG sample. 
This study also revised the redshift of LOCK850.01 from 
z - 2.148 to z - 3.38. We adopt this new redshift because 
previous estimates were based on the spectroscopic follow-up 
of a galaxy ~ 3" away from the radio counterpart of this submm 
source. 



The resulting sample of 18 SMGs with robust redshift 
estimates was cross-matched with our MIPS-PACS-SPIRE 
multi-wavelength catalogue. Fifteen are detected in at least 
one of the PACS/SPIRE bands. Tables |7] and [8] present the 
multi-wavelength properties of this subsample. 



3.4. COSMOS 

In the COSMOS field we use the multi-wavelength identifica- 
tion of LABOCA and MAMBO sou r ces ca rried out by Aravena 
et al. (in prep.) and i Bertoldi et al.l (|2007|) . respectively. From 
the Aravena et al. sample we only keep sources with radio 
identifications. This limits our sample to 46 SMGs out of the 
163 LABOCA sources. In the Bertoldi et al. sample there are 
27 MAMBO sources with robust radio identifications. Among 
those sources, nine are already included in the Aravena et al. 
sample. For those sources we keep the radio identification 
obtained by Aravena et al. because it is based on the latest 
version of the deep COSMOS radio catalogue. 

We cros s-match this sam ple of 64 SMGs with the AzTEC 
catalogue of IScott eFd ] (|2008), which has no multi-wavelength 
identifications. AzTEC sources with no LABOCA or MAMBO 
counterparts but wit h Submillimeter Array (SMA) follow-up 
(I Younger et al.ll2007i l2009i) are included in our sample (i.e., 5 



sources). 

Capak et al. (in prep.) provide redshift follow-up for some 
of these 69 SMGs with robust multi-wavelength identifications. 
So far this spectroscopic follow-up programme has obtained 
redshift estimates for 15 of these SMGs. 

These 15 SMGs with robust redshift estimates are 
cross-matched with our MIPS-PACS-SPIRE multi-wavelength 
catalogue yielding 11 SMGs detected in at least one of the 
PACS/SPIRE bands. Tables |9] and [10] present the multi- 
wavelength properties of this subsample. 



3.5. Cluster fields 

We gather from the literature a sample of well-known lensed 
SMGs with both spectroscopic redshifts and lensing magnifica- 
tion e stimates. In the A22 18 fi eld, our SMG sample i s assem bled 
from iKneib et all (|2004 and iKnudsen et al] (120061 12008 ) and 
contains six lensed sources. Among these six lensed sources, 
three correspond to the same lensed galaxy (SMMJ16359-I-6612; 
Kneib et al. '"200?). In A1835, submm observations are taken 
from Ivison etal. (200( 1). The re dshift of SMMJ1401 1+0252 
is also taken from Iviso n et alj (120001) , while the redshift 
estimate of SMMJ14009H-0252 is from IWeiB et all (l2009al) . 
In MS0451 and A22 19, submm observations are taken from 
Chapma n et al.l(l2002l) . Each field contains only one lensed SMG 
with both spectroscopic redshifts and lensing magnification 
estimates, namely, SMMJ 16403 +4644 and SMMJ04554+0301 
(Rigbv et al.t,2008; .Borvs et al...2004. respectively). In MS 1054, 
we use submm observati ons and redshift information provided 
in 'Knud sen et aO (200^', SMMJ10570-0336). For A1689, we 
use s ubmm observat ions and lensing magnification estimates 
from iKnudsen et all (l2008l SMMJl 31 15 - 1208) while red- 
shift informations are from togbv et"aD (l2008h . Finally in 
CL0024, A2390 an d A370 submm observations are taken from 
Small et al.' ("2002). The re dshift of SMMJ00266+1708 comes 
from Vahante et all (l2007h . the redshift of S MMJ02399-0136 
comes from llvison et al.l (11998) (see also iLutz et al] 2005h 
and the redshift of SMMJ02399-0134 comes from [S mail et ajj 
(2002). For SMMJ21536+1742 we use lBarger etal] ( 1999) (K3 
counterpart; Fraver et a l. 2004). 

All but one of these sixteen lensed SMGs have been 
detected in at least one of the PACS/SPIRE bands. Because 
these galaxies are magnified, their mid-to-far infrared fluxes 
are de-magnified prior to further analysis using magnification 
factors from the Uterature. Tables [TTI and fT2] present our lensed 
SMG sample. 



The infrared luminosities of our lensed SMGs strongly de- 
pend on their magnification factors. These factors are estimated 
from complex lens models, constrained by the many lensed 
features seen in these clusters. We adopt a characteristic error 
of 20% on their luminosities to account for uncertainties in the 
lens models. 



3.6. SMGs with multiple counterparts 

Our SMG sample contains 62 sources detected by PACS/SPIRE 
and with secure spectroscopic redshift estimates. Among 
these 62 SMGs, eleven have multiple optic al/radio/MIPS 
counterparts. Six of them (GN04, GN07, GN19, GN39, 
AzTECJ 100008 +024008 and MAMBOll) are treated as one 
single system because they are assumed to be interacting galax- 
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ies. The optical counterparts of GN19 and GN 39 are spectro- 
scopically confirmed to lie at the same redshift ("Chapm an et alJ 
[2005; Swinbank et al. 2004) and the optical counterparts of 
GN04 and GN07 exhibit IRAC photometry consistent with both 
optical sources being at the same redshift. The optical coun- 
terpart of MAMBOll without any spectroscopic redshift es- 
timate (MAMBOllW) has a photometric redshift supporting 
the assumption of an interacting system (Bertoldi et al. 2007J. 
AzTECJ 100008 -1-024008 has two SMA c ounterparts within th e 
submm beam with consistent redshifts dYounger et al.l l2009h . 
Because these multiple counterparts are thought to be part of an 
interacting system, to derive the dust properties of these galaxies 
we sum the mid-infrared, far-infrared and radio flux densities of 
their optical/radio/MlPS counterparts. 

For four SMGs we have a spectroscopic follow up for only 
one of their multiple MlPS/radio counterparts, LOCK850.03, 
LOCK850.04, LOCK850.15 and LESS 10, namely. Thus we can- 
not assess whether these galaxies are interacting systems. We as- 
sume that only the source with a redshift estimate significantly 
contributes to the submm and far-infrared flux-densities. This as- 
sumption is supported by the fact that the MlPS-24 fim and radio 
flux densities of these sources agree with the infrared luminosi- 
ties derived from their far-infrared/sub mm flux densities. The 
inclusion or exclusion of these four sources would not change 
the conclusions of our study. 

LOCK850.41 has two robust radio counterparts coincid- 
ing with two MlPS-24yum sources. Spectroscopic follow-up 
of these counterparts shows that they do not correspond to 
an interacting system, one galax y is situated at z = 0.689 
("Menendez-Delmestre et al.l 12009 "') and the other at z = 0.974 
(Coppin et al. 2010). IRS observations show that while the low 
redshift galaxy exhibits strong PAH signatures, the galaxy sit- 
uated at z = 0.974 has a continuum-dominated mid-infrared 
spectrum with no visible PAH features, consistent with an AGN 
classification. This suggests that the high-redshift galaxy has 
very low ongoing star-formation, incompatible with bright far- 
infrared and submm emission. However, because this assump- 
tion is still highly uncertain, we decide to remove this source 
from our final sample. 

3. 7. Stellar mass estimates 

Due to the significant obscuration at rest-frame optical wave- 
lengths, and to the possible presence of a rest-frame nea r-lR 
continuum excess in numerous SMGs (iHainline et al.ll20Tlh . the 
determination of the stellar mass of SMGs is still highly debated. 
For example, different assumptions about the star-formation his- 
tory or about the contribution of an AGN to the rest-frame near- 
IR continuum excess could lead to systematic variations in the 
med ian stellar mass estimates of SMGs of more than a fa ctor 2 
(see iHainUne et all 120111: iMichalowski et al.ll2010l 1201 i V Due 
to all these different methods and assumptions, it was impos- 
sible to find stellar masses homogeneously derived for all our 
SMGs in the literature. Therefore, we decided to infer the stellar 
masses of our SMGs using a single method. We would like to 
stress that resolving the problem of the stellar mass estimates of 
SMGs is beyond the scope of this paper The absolute values of 
our estimates might not be fully reliable, but the fact that we are 
using a homogeneous method and assumptions over our sample 
should provide a good tool to study relative variations in stellar 
mass. Lensed SMGs are not considered in that study because of 
the difficulty to obtain coherent optical-to-near infrared data for 
these galaxies, making any stellar mass estimates very uncertain. 
Optical-to-near-infrared photometry was obtained using the 



radio or optical positions of our SMGs. In GOODS-N and 
COSMOS, we used the multi-wavelength catalogue built by 
the PEP consortium and presented in Berta et al. (2010, 20lib. 
In GOODS-S, we used the MUSIC catalogue (S antini et d] 
2009) and t he optical-to-near-infr ared photometry of SMGs 
presented in IWardlow et al.l (1201 lb . In the LH field, we used 
the optical-to-near-infrared photometry of SMGs presented in 
Dye et al. (2008) and Copp in et al. (2010). Stellar masses were 
then calculated by fitting the multi-wavelength photometry to 
iBruzual & Charlon (l2003h te mplates through a m inimization, 
usin g the method describ ed in lFontana et al.l (l2004t) and updated 
as in lSantini et al.l (12009*). We looked at all fits and rejected those 
sources with problematic fits. Among the 46 SMGs considered 
in this study (all our blank field SMGs), 39 SMGs have good 
optical-to-near-infrared SED fits. The stellar masses of these 39 
SMGs are provided in Table [13] In the LH field, we find that 
our stellar mass estimates are in perfect agreement with results 
from Hainline et al..(.201 !■) . The agreement between our findings 
is encour aging because CO observations and dynamic mass ar- 
guments (lEngel et al.ll2010h favour these lower stellar mass esti- 
mates, more consistent with the finding s of Hainline et al. (201i]) 
than those of Michalo wski et al ] (12010 ). The median log(M,) of 
10.86 for our sample is also fully consistent with l og(M,) ~ 1 1.0 
for SMGs estimated from the SMG halo mass of Hickox et al. 
(20121), using the conversion to stellar mass by iMoster et al. 
(20m . 

3.8. Final sample and selection biases 

Our final SMG sample contains 61 sources detected by 
PACS/SPIRE and with secure spectroscopic redshift estimates. 
Because this sample requires MIPS detections, PACS or SPIRE 
detections and robust redshift estimates, it is aff'ected by several 
selection biases. Previous studies have already discussed the 
biases introduced by (sub)mm o bservations and/or SPIRE-like 
(i.e., BLAST) observations (e g., iCasev et alj2009tlChapin et al.1 
1201 iblSvmeonidis et alJl20ll but none of them have examined 
our peculiar selection function. In this section we list all our 
selection biases and try to estimate how representative our 
sample is of the SMG population and more generally of the 
high-redshift star-forming galaxy population. Here, we only 
focus on the blank field SMG population because lensed SMGs 
are affected by more complex selection function depending on 
their positions with respect to the foreground lenses. 

Because (sub)mm and far-infrared surveys observe 
the thermal emission of dust they are limited, at a given redshift, 
in the range of infrared luminosities and dust temperatures 
probed. In order to quantify these selection biases we studied 
the Tdust - Lib. parameter space reachable with our far-infrared, 
submm and radio observations. For that purpose we took a 
model describing the far-infrared SED of SMGs (a power-law 
temperature distribution parameterized with T^, i.e., the temper- 
ature of the coldest dust component of the model, see section 
14.21 ) and estimated for each point of the - Lir parameter 
space its detectability by the PACS (lOOyum or 160jum but 
mainly by the 160/im band), SPIRE (250 fim, 350 yum or 
500 yum but mainly by the 250 //m band) and SCUBA (850 yum) 
instruments. Then, in order to com pare these estirn ates with the 
local Tdust - Lib relation derived bv lChapman et alj (|2003) using 
a single temperature optically thin modified blackbody model, 
we simply converted into Tdust with - 0.6 x T^asi + 3 K 
(see section 14.2.21 and Fig. |5j. This study cannot be directly 
performed using a single temperature optically thin modified 



Magnelli et al.: Far-infrared properties of SMGs 



7 




Fig. 1. {Left) Selection limits introduced in the Tjust - Lir parameter space by single-wavelength detection techniques. Continuous, dashed, 
dotted-dashed, triple-dotted-dashed and dotted lines show the lower limits on ZjjR introduced by the submm, PACS, SPIRE, MIPS-24 jt/m and radio 
observations, respectively, at z = 1.5 (thin blue lines) and at z = 2.5 (thick red lines). The parameter space reachable by a given single-wavelength 
detection technique corresponds to the area situated to the right of the lines. As an example, the red arrows show the para meter space probed at 
z ~ 2.5 by our GOODS-N SMG sample. The shaded area shows the loc al Tdust ~ ^ir r elation found bv lChapman et al.l ( l2003i) . linearly extrapolated 
to 10'^ Lq. The striped area presents results for SMGs extrapolated bv lChapman et al. (2005) from radio and submm data. {Right) The hatched 
histogram shows the redshift distribution of our PACS/SPIRE detected SMG sample. The empty histogram shows the redshift distribution of its 
parent sample, i.e., SMGs with robust redshift estimates obtained through secure multi-wavelength identifications. 



blackbody function because that model cannot reproduce the 
PACS 100//m measurements sometimes dominated by warmer 
or transiently heated dust components (see Section 14.1b . For 
the radio detectability we used the local far-infrared/radio 
correlatioiQ (Helou et all 119881: lYun et alj 1200 Ih and f or the 
MIPS-24 um detectability we used the ICharv & Elbad (1200 lb 
template^ In this exercise we used the typical 3cr limits of 
GOODS-N observations, i.e., 20juJy, 3mJy, 6mJy, lOmJy, 
12mJy, 12mJy, 3mJy and 15//Jy at 24 /vm, lOOjum, 160/vm, 
250 //m, 350 jum, 500 jum, 850 yum and 1.4 GHz, respectively. 
The left panel of Fig. [1] shows the selection limits observed in 
the GOODS-N field. To obtain the selection functions of the 
other fields, one would simply shift the lines of Fig. [T] towards 
higher infrared luminosities according to the depth of the 
observations with respect to the GOODS-N field (see Table|2]l. 

The first selection bias introduced in our SMG sample 
come from the (sub)mm detections. This selection bias is almost 
redshift independent, but selects, at a given infrared luminosity, 
only galaxies with cooler dust. The bias decreases at high 
infrared luminosities where submm observations probe a large 
range in dust temperature. In fields where (sub)mm observations 



' In Section 15.11 we find that the parameterization of the far- 
infrared/radio correlation, {q), is slightly lower in our SMG sample than 
in the local universe, {q) = 2.0 versus (q) = 2.34. However, here, we 
prefer to use the local value of (q) because our sample cannot be used to 
fully constrain this parameter. This is a conservative approach because 
using a lower value of {q} one would decrease the selection bias intro- 
duced by radio observations, i.e., radio observations could reach lower 
infrared luminosities at a given redshift. 

In Section 15.11 we find that the infrared luminosities estimated 
from the MIPS-24 ;/m fluxes densities and the Charv & Elbaz (2001) 
library are ove r estim ated. Therefore, in this exercise, the use of the 
ICharv & ElbazI ( 1200 D library is a conservative approach because at 
high-redshift the MIPS-24 fim observations could reach even lower in- 
frared luminosities. 



are shallower than in GOODS-IS0, these selection functions 
shift towards higher infrared luminosities. Nevertheless, shallow 
(sub)mm observations would still probe, at high infrared lumi- 
nosities, a large dynamic range in dust temperature. Therefore, 
assum ing that the local Tdu s t - Lir relation holds at high redshift 
(e.g. iHwang et all 120101: IChapin et all 120111: iMarsden et all 
1201 1*), and extrapolating it towards higher infrared luminosities, 
we can assume that at high luminosities (Lir > lO'^'^'L©), 
SMGs are a representative sample of the underlying star- 
forming galaxy population. 

The second selection bias affecting our sample comes from 
the necessity of having robust redshift estimates. This require- 
ment translates into accurate positions and multi-wavelength 
identifications mainly obtained via radio observations (among 
the 69 SMGs with redshift estimates in our blank fields, 59 have 
been identified using radio observations while only 5 have been 
identified using MIPS-24 jjm observations and 5 using SMA ob- 
servations). Radio observations probe the synchrotron emission 
of galaxies and suffer from positive A;-corrections, independent 
of the dust temperature. This biases our sample towards higher 
infrared luminosities as the redshift increases (see dotted lines 
in the left panel of Fig.[T]i- The redshift estimates of these 
radio sources, obtained mainly through optical spectroscopy, 
introduce additional selection biases. For example, just for 
feasibility of the optical spectroscopy in a reasonable amount 
of time and/or success of detection, SMGs with spectroscopic 
redshifts mi ght be biased towards optically-bright SMGs (see 
e.g.. Chapman et al.ll200 5h and are also likely to have a higher 
incidence of strong emission lines than typical SMGs. In 
addition, spectroscopic follow up of SMGs might also miss 
some objects at 1.2 < z < 1.8 (namely the "redshift desert"), 
due to the lack of strong emission lines in the rest-frame 
wavelength range o bserved by ground-b ased spectroscopic 
instruments (see e.g.. IChapman et al]l2005h . All these selection 



' One can convert the MAMBO or AzTEC flux density limits into 
its corresponding SCUBA-850 //m flux density limit using the Raleigh- 
Jeans approximation. 
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biases are very difficult to quantify because they depend on the 
follow-up strategy used. Here, using a Kolmogorov-Smirnov 
(KS) analysis, we simply verify that the radio and submm flux 
density distribution of SMGs with spectroscopic redshift is 
consistent with that of its parent sample, i.e., SMGs with radio 
counterparts. This suggests that the spectroscopic follow-up of 
radio-identified SMGs does not introduce strong biases towards 
any particular infrared luminosity or dust temperature. On 
the contrary, we find that the distribution of submm to radio 
flux ratio of the SMGs with spectroscopic redshift is different 
than that of its parents sample (only 30% of chance for being 
drawn from the same distribution). Because the submm to 
radio flux ratio has been used as a reds hift indicator by many 
early works (e.g., ICariUi & Yuril Il999t IChapman et alj|2005h . 
we conclude that spectroscopic follow-up of SMGs might 
be slightly biased towards low redshift galaxies. However, in 
terms of luminosities and dust temperatures, we assume that at 
high infrared luminosities (Lir > 10'^'^ L©), SMGs with robust 
spectroscopic redshift estimates are still a good representation 
of the underlying SMG population and therefore of the entire 
high luminosity star-forming galaxy population. At low infrared 
luminosities, however, SMGs with redshift estimates represent 
a subsample of SMGs biased towards lower redshift galaxies, 
essentially because of the need for a radio-based identification. 

Our final SMG sample is also affected by the MIPS-PACS- 
SPIRE detection requirement. The MIPS-24yum requirement 
should not significantly influence our sample because it corre- 
sponds, up to z ~ 3-4 and in all our fields, to selection limits 
several times lower in term of infrared luminosities than those 
introduced by radio observations (see triple-dotted-dashed line 
in the left panel of Fig. [B- On the contrary, the PACS/SPIRE 
requirement affect our sample and is redshift dependent. PACS 
observations, which suffer from positive A:-corrections, are 
slightly biased towards galaxies with hotter dust while SPIRE 
observations are biased towards cooler dust. The SPIRE selec- 
tion bias is also redshift dependent because SPIRE detections 
are mainly obtained in the 250 /im band which suffers from 
positive ^-corrections as it reaches the peak of the far-infrared 
SED of galaxies at z~1.5. In GOODS-N, the selection bias 
due to the PACS/SPIRE observations is almost equivalent 
to that introduced by the combination of submm and radio 
observations. In other fields, the PACS/SPIRE requirement is 
even less constraining because the SPIRE observations are as 
deep as in GOODS-N while radio and (sub)mm observations 
are shallower This is reflected by the fact that the PACS/SPIRE 
detection rate of SMGs with robust spectroscopic redshift 
estimates is very h i gh, an d much higher than that observed by 
iDannerbauer et all (l2010h for the entire SMG population, i.e., 
73% versus 39%. 

In summary, our final SMG sample should provide a good 
representation of the high infrared luminosity (Lir > 10'^'^ Lq) 
SMG population and more generally, of the entire high infrared 
luminosity galaxy population. On the other hand, as we go to 
lower infrared luminosities (Lir < 10'^'^ Lq), our final SMG 
sample is biased towards low redshift galaxies with cold dust. 
Most of these biases are not inherent to our PACS/SPIRE 
SMG subsample but are intrinsic to any SMG sample requiring 
robust spectroscopic follow-up aided by secure radio/MIPS 
multi-wavelength identifications. 

The right panel of Fig.[T] presents the redshift distribu- 
tion of our PACS/SPIRE-detected SMG sample. This redshift 
distribution is consistent with that of the entire SMG sample 



with robust redshift estimates. The median redshift of our 
PACS-SPIRE detected SMG sample is z = 2.4 and is consistent 
with th e median redshift of the entire SMG population, i.e., 
z ~ 2.3 (IChapman et al.ll2005h . 

4. SED analysis 

In this section we describe the models used to infer the dust 
properties of the SMGs. Scientific conclusions drawn from 
these properties are discussed in Section|5j 



4. 1 . Single modified blackbody model 

In order to infer the dust temperatures and infrared luminosities 
of our galaxies we fitted their far-infrared and (sub)mm flux 
densities with a single temperature modified blackbody model. 
This model provides a very simple description of the far-infrared 
SED of a galaxy, because it assumes that the emission-weighted 
sum of all the dust components could be reasonably well 
fitted by only one blackbody function at a given temperature. 
Despite its simplicity and the fact that it is known that this 
model cannot fully reprod uce the Wien side of the far-infrared 
SED of galaxies ( e.g., .Blain et al.ll2003t iMagnelU et aP 120101: 
iHwang et alJl2010l) . we adopted this model for two reasons: (i) 
studies of the Infrared Astronomical Satellite (IRAS) galaxies 
have demonstrated that it still provides an accurate diagnostic 
of the typical heating conditions in the in terstellar medium of 
big grains in thermal equilibrium jDesert et al., 1990) : and (ii) it 
allows direct comparison with most of the pre-Herschel studies. 
The far-infrared flux densities of our galaxies were thus fitted, 
in the optically thin approximation, with a single modified 
blackbody function : 

SyCC , (1) 

exp(/!v/^rdust) - 1 

where 5v is the flux density, /3 is the dust emissivity spectral 
index and Tdust is the dust temperature. This single temperature 
modified blackbody model cannot reproduce the full rest-frame 
8-to-1000;um SED over which the total infrared luminosities 
(LiR[8-1000//m]) are classically defined. A significant amount 
of energy emitted at relatively short rest-frame wavelengths 
(i.e., where the backbody function drop sharply) would thus 
be missed by a simple integration of the blackbody function 
over the rest-frame 8-to-1000//m wavelengths. Therefore, 
the total infrared luminosities of our galaxies were inferred 
using the far-infrared luminosity definition (LFiR[40-120^(m]) 
given by Helou et al. (1988) and a bolometric- correction term. 
This bolometric-correction is equal to 1.91 ( Dale et al.ll200U 
LiR - 1.91 X Lfir) but introduces uncertainties in our estimates 
becau se it varies (±30 %) with the intrinsic shape of the galaxy 
SED (iDale et al.ll200l . 



4.1.1. Constraints on /? 

The exact value of the dust emissivity spectral index /3 is still 
debated. Laboratory experiments as well as observations in di- 
verse Galactic environments suggest a broad range of values 
for /3 (Dunne & Eales 2001; Dupac et al. 2003, and references 
therein). The value of /3 seems to depend on the chemical com- 
position, the temperature and the size of the dust grains. Despite 



Magnelli et al.: Far-infrared properties of SMGs 



9 



its variability on Galactic scales, extragalactic constraints on yS 
converge to a narrow range of values (1.5 < /} < 2.0). In par- 
ticular, ODunne & Bales (2001) found a constant dust emissivity 
spectral index of ~ 2 using a sample of galaxies probing a 
broad range of infrared luminosities. Based on this latter conclu- 
sion, we assume that /? could be considered as universal over the 
SMG population. 

Assuming (3 to be universal, we can constrain its value glob- 
ally using our sample of 61 SMGs. To perform this global fit 
we gridded the p parameter space [0.1-3.0] with steps of 0.05. 
Then, for each value of y6, we performed a minimization for 
each galaxy, varying Tdust and the blackbody normalization. The 
value at a given fi is then defined as the sum of the;\f^ value of 
all galaxies (i.e., x^. - 2 xl^i)- Our;^^ minimization was done 
using a standard Levenberg-Marquardt method. 

We apply this global fit to three different wavelength cover- 
ages. First, we fit the full wavelength coverage provided by the 
Herschel and (sub)mm observations (i.e., from the PACS 70jum 
to the (sub)mm wavelength); second, we exclude from the fits 
the PACS 70 and 100 yum data points; and third, we exclude from 
the fits the PACS 70, 100 and 160yL(m data points. For these 
three different wavelength coverages the best fit is obtained at 
j8 = 0.6 + 0.2, p = 1.2 + 0.2and^ = 1.7 + 0.3, respectively (using 
the 95% confidence level, i.e., Ax^^x^ia ^'^^^ ^^'^^ these 

errors stand for the mean values, rather than for the standard de- 
viation of the population). Fits of the full wavelength coverage 
systematically lead to significantly larger values than for the 

other cases (i.e., A'gai~18 for Ndof~4). On the contrary, fits ex- 
cluding the PACS 70 and 100 yum data points or the PACS 70, 
100 //m and 160 fim data points lead in both cases to low xl^i 

values, i.e., v/ithxl^^ ~ 6 for Ndof ~ 3 and;^^^^; ~ 4 for Ndof ~ 2, 
respectively. 

The large values observed when we try to reproduce the 
full wavelength coverage provided by the Herschel and (sub)mm 
observations perfectly illustrate the limits of a single temper- 
ature model. Such a simple model cannot fully d escribe the 
Wien side of the far-infrared SED of galaxies (e.g..'Bla in et alj 
[2003; MagnelH et al. 2010; Hwang et al. 2010). The PACS 70 
and 100 /im flux densities are likely dominated by a warmer or 
transiently heated dust component. Consequently, the PACS 70 
and 100 fj.m data points have to be excluded from the fitting pro- 
cedure. A precise description of the far-infrared SEDs of galax- 
ies requires a more complex model which includes multiple dust 
components (see Section l4~2l i. 

The increase of /3 when excluding short-wavelength 
measurements fro m the fits agrees with the conclusions of 
IShettv et al.l (j2009) studying galactic dense cores : constraints on 
/3 are highly sensitive to the wavelength coverage used in the fits 
as well as to the noise properties of the observations. Although 
interesting, our constraints on /3 should thus be used with cau- 
tion. 



4.1 .2. Fitting tine full SMG sample 

In the following, we decide to fix the dust emissivity spectral 
index /3 to its standard value of 1.5. This choice is driven by two 
reasons. First, this value is fully compatible with our findings 
(i.e., 1.2 < /3 < 1.7) and second, it allows direct comparison with 
all pie-Herschel studies. We also decide to exclude from our 
fits the PACS 70 and 100/vm data points because they are likely 
dominated by a warmer or transiently heated dust component. 
The PACS 160yum data points are kept because their exclusion 



does not significantly improve our fits while their inclusion 
allows better constraints of the dust temperature estimates (see 
Fig. 13. 

Figure IATI presents results of this fitting procedure to each 
individual SED, while Table [13] gives the inferred dust temper- 
atures and infrared luminosities. Uncertainties are estimated 
using the distribution of Tdust and Lir values that correspond to 
models with < ■ + 1 . 

^ ^min 

We observe in Fig. lA.ll that a single dust temperature model 
provides a reasonable fit to the 160//m-to-mm data points 
(with x^ ~ 7 for Ndof ~ 3). Figure [ATI also shows the limits 
of this model at short wavelengths and why we excluded from 
our fits the PACS 70jum and 100 /vm data points. The modified 
blackbody functions drop quickly at short wavelengths and 
cannot reproduce the PACS 70 //m and lOOyum data points of 
most of our SMGs. 

For some of the SMGs we do not have both PACS 
and SPIRE detections. For those galaxies, we can expect the 
inferred dust temperatures and infrared luminosities to be more 
uncertain, and potentially biased because PACS and SPIRE 
measurements probe different parts of the blackbody emission 
of the dust (Wien and Rayleigh-Jeans side, respectively). To 
assess this issue, we compared the dust temperatures and 
infrared luminosities inferred using the combination of PACS 
and submm observations, or SPIRE and submm observations, to 
the reference values inferred using the continuous wavelength 
coverage provided by the combination of PACS, SPIRE and 
submm observations. This analysis is based on 50 SMGs 
detected by both PACS and SPIRE. Results are shown in Fig.|2] 

For most of our sources the dust temperatures and in- 
frared luminosities estimated from the combination of PACS 
(or SPIRE) and submm observations are in good agreement 



PACS 



with our reference values, i.e., a-\Tf'^\ - T, , 

' ' dust dust 



((r[T^ 



Ref 



- T 



SPIRE 



] = 2.5 K 
] = 3.1K) and cr[L!l'^VLr^'^^] = 0.08 dex 



(cr[L|^*7L^P''*^] = 0.1 dex). However, the dust temperatures 
inferred using SPIRE and submm observations are slightly 
underestimated at high dust temperature (Tdust > 35 K). At 
these temperatures, the SPIRE observations start to be affected 
by the shift of the far-infrared SED peak towards rest-frame 
wavelengths barely probed by the SPIRE 250 yum passband. 
This effect slightly biases these estimates. 

There are only a few sources with large uncertain- 
ties (i.e., Ar>8K or Alog(LiR)>0.3 , COSLA127R1I, 
AzTECJ100019+0232, SMMJ105238+5716, GN26, and 
SMMJ163541+6611). Examining the SED fits of these galax- 
ies, we find that all of them exhibit large x^ (i-^-, Z 15 for 
Ndof~3) when combining their PACS, SPIRE and submm 
observations. These large x^ values seem to be explained by 
one or two inconsistent flux densities in their SEDs. These in- 
consistent data points do not correspond to a specific rest-frame 
wavelength but randomly affect the PACS, SPIRE or the ground 
based data points. Thus they are unlikely due to st rong e mission 
lines (like the [C II] emission line. Small et al. 1201 ih which 
are not included in our simple modified blackbody model. We 
conclude that the observed discrepancies are not directly due to 
our simple modified blackbody model but to some outlying flux 
densities, as expected when working close to the non-Gaussian 
confusion limit which can create significant outliers. 

Finally, one can expect the accuracy of the estimates 
inferred from the combination of PACS (or SPIRE) and submm 
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Fig. 2. (Le/?) Dust temperatures inferred from the combination of PACS only (or SPIRE only) together with submm observations, compared 
with the reference values inferred using PACS, SPIRE and submm observations. These comparisons are for a single dust temperature modified 
blackbody model. Blue squares represent SMGs situated in blank fields while green diamonds represent lensed-SMGs. {Right) Same comparison 
but for the inferred infrared luminosities. The dust temperatures and infrared luminosities of galaxies can be reasonably inferred from their 
PACS+submm or their SPIRE -l-submm observations alone using a single temperature modified blackbody model. 



observations to vary as function of the redshift: high(low) 
redshift galaxies with PACS (SPIRE) only measurements 
could have inaccurate dust temperature estimates because their 
far-infrared SED peak shifts outside the PACS (SPIRE) bands. 
However, we find no significant evolution of Ar or Alog(LiR) 
with the redshift. At low redshift, the shift of the far-infrared 
SED peak towards shorter wavelengths is counterbalanced by 
the fact that at these redshifts, galaxies exhibit relatively low 
infrared luminosities and dust temperatures, shifting back their 
far-infrared SED peak towards the SPIRE bands. Likewise, at 
high redshift, SMGs exhibit higher infrared luminosities and 
dust temperatures, shifting back their far-infrared SED peak 
towards the PACS bands. 

We conclude that the dust temperatures and infrared 
luminosities of galaxies can be reasonably inferred from their 
PACS-Hsubmm or their SPIRE -Hsubmm observations alone using 
a single temperature modified blackbody model. This may be 
important for survey regions covered at sufficient depth with 
one of these instruments only. 

4.2. Power-law temperature distribution 

Although a single-temperature model gives a good description 
of the far-infrared peak and Rayleigh-Jeans side of the SED 
of SMGs, it fails to reproduce short wavelength observations 
(e.g., the PACS 70 and 100 fim passbands) which are affected by 
warmer or transiently heated dust components. Consequently, 



the total infrared luminosity of SMGs (i.e., LiR[8-1000/im]) 
has to be extrapolated from LFiR[40-120yL(m] and short wave- 
length observations have to be excluded from the fit. In order 
to reproduce these short wavelength observations we need to 
use a more complex model, taking into account warmer dust 
components. 

T o de scribe the dus t emission of galaxies, iDale et all 
(1200 Ih and iDale & Heloul (2002) assumed that they are the 
superposition of regions heated by different radiation fields. In 
that framework, they assumed that the dust mass submitted to 
a radiation field U is given by dM^ast/dU oc U^". Then using 
simple assumptions they showed that a ~ 2.5 is appropriate 
for a diffuse medium whil e g ~ 1 descr ibes^ a dense medium. 
Following the same idea, iKovacs et al.l ilOlO) described the 
SEDs of galaxies by a power-law distribution of temperature 
components (dM^ust/dT oc T^'>') with a low-temperature cutoff 
Tc. Under the assumption that the dust is only heated by radia- 
tion (and not by non-radiative processes like shocks), the main 
parameters of this model and that of Dale & Helou are linked 
hy y X 4 + a + ySeff (where ySetf is the dust emissivity spectral 
index observed near the peak of the far-infrared emission). This 
model can accurately describe the mid-to-far-infrared SEDs of 
local starbursts (Kovacs et al. 2010) and is convenient for our 
purposes as it is parameterized in dust temperature rather than 
radiation field. Consequently, while other models could have 
been used (e.g., iDale & Heloul 120021: iDraine & Lil l2007h . we 
adopted this prescription as a natural extension of our single 
dust temperature model. 
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The parameterization of this power-law temperature 
distribution model is fully described in Kovacs et al. (2010), 
and briefly summarized here. In particular we do not give 
the analytical derivation of the infrared luminosity because 
here we derive this quantity using a simple discrete numerical 
integration. 

Expressed in observable parameter space, the emission 
from a single modified blackbody emission, not in the optically 
thin approximation, is given by 



(2) 



where By is the Planck function, Tobs is the observed-frame tem- 
perature (i.e., Tobs = T/{1 + z)), T is the optical depth, dQ. is the 
solid angle subtended by the galaxy and m is a magnification cor- 
rection for lensed galaxies (= 1 in all other cases). In the model 
proposed by Kovacs et al. (2010), the optical depth is expressed 
as a function of the dust mass (Mdust) and the projected source 
diameter (R), together with the usual power-law frequency de- 
pendence for the emissivity of dust. 



T(Vr) = Ko 



ttR^ 



(3) 



where t is expressed in the rest-frame (v,- - Vob.s(l + z)) 
and K() is the photon cross-section to mass ratio of parti- 
cles at t he reference f requen cy vq. To allow direct compari- 
son with l^vicsietaL j 201^2_we adopted /cgso - 0.15 m^ kg ' 
at vo = c/850//m "(•Dunn e et alJ l2003h . even though the ex- 
act value of this parameter is still und er active discussion 
(e.g., Eil debrand 1983; Ki-ugel et al. 199d: ISodroski et al.ll997t 
iJames et al.,.2002) . Using this formalism a power-law tempera- 
ture distribution model can be expressed as. 



where is the low-temperature cutoff of the model. 



(4) 



4.2.1 . Constraints on j3, y and R 

The power-law temperature distribution model has five free 
parameters, Tc, Mdust, J and 7?. It can only be constrained from 
observations that probe the full far-infrared SEDs of galaxies, 
i.e., probing the Wien-side, the peak and Rayleigh-Jeans-side 
of these SEDs. Such broad spectral coverage can only be 
obtained through the combination of PACS, SPIRE, submm and 
millimeter observations and thus can only be applied to a small 
fraction of our SMG sample. Therefore here we investigate the 
possibility that some of those parameters are universal over the 
full SMG population. 



As akeady mentioned in Section 14.1.11 considering 
that the exact value of the dust emissivity spectral index fi is 
still debated, one can assume this value to be universal over the 
SMG population. 

iKovacs et alJ (l2OI0l) found little variation of y in the 
local star-forming galaxy population. Based on this finding 
they assumed a constant value of j for high-redshift luminous 
starbursts and obtained a good fit to their SEDs. Therefore, 
in the following, we consider y as universal over the SMG 
population. 

Finally, we considered the projected radius of the emitting 



region, R, as universal over the SMG population. This consid- 
eration is perhaps questionable because in high-redshift star- 
forming galaxies, the diameter of the region forming stars spans 
a wide range o f values from 1 to 10 kpc (Chapman et al. 2004t 
'M uxlow et al.ll2005l: iTacconi et al.ll2006> ,2008; .Biggs & IvisonI 



200 81: ICasev et al.1 120091: llono et alj |20()9|; iLehnert et'aiT 
CarilH et al.' '2010: Swinba nk et al.1 120IOI: iTacconi et al.1 



2009; 



2O10; 



Younger etal. 2010). However, in the power-law temperature 
distribution model the variation of R does not strongly affect 
the estimates of Ljr (< 5%) but only affects the physical 
interpretation that one can draw from the absolute value of T^: 
smaller values of R imply higher values for T^. In any case, the 
study of the relative variation of from one galaxy to the other 
is not qualitatively affected by the exact value of R. 

Assuming these three parameters to be universal, we con- 
strained them globally using a subsample of 19 SMGs detected 
in all PACS and SPIRE passbands and with at least one detection 
longward of 1 mm (needed to obtain good constraints on the dust 
emissivity /?). To perform this global fit we first gridded the y 
and R parameter space using ranges of [1.0-2.5], [6.5-9.0] and 
[0.5 kpc-9.0kpc] and steps of 0.05, 0.1 and 0.25, respectively; 
then, for each node of this grid, we performed Ilx^ minimization 
for each galaxy, varying and Mdust- The value of the 
node is then defined as the sum of the value of all galaxies 
{i.Q., x\aAtt ~ 2 ^gai)- Our;^'^ minimization was done using a 
standard Levenberg-Marquardt method. Figure [3] presents the 
confidence levels obtained for fi, y and R. Confidence levels are 
computed using Kx'^^x]^,,^ + [2.3,6.0, 11.6] for the 68%, 95% 
and 99% confidence level, respectively. The best fit is obtained 
at yS = 2.0 + 0.2, y = 7.3 + 0.3 and = 3 + 1 kpc (using the 
95% confidence level; note that these errors stand for the mean 
values, rather than for the standard deviation of the population), 
and corresponds to;\fg^j~8 for Ndof ~5. These values confirm 
that our model provides a good description of the far-infrared 
SEDs of SMGs even if three parameters are considered common 
to all galaxies. 

In Fig. [3] we observe only small degeneracies between yS, 
y and R, e.g., an increase of ji could be compensated, in terms 
of x^ minimization, by an increase of R. The wide wavelength 
coverage provided by our data allows us to reasonably constrain 
our model. Constraints on fi, y and R are also in line with 
the physical expectations and with independent estimates. A 
dust emissivity spectral index j6 of 2.0 + 0.2 is in agreement 
with conclusions based on local L IRG/ULIRG dDunne & EalesI 
1200 ll: IChakrabarti & McKeell2008 ) . The dust emissivity spectral 
index found using our power-law temperature distribution 
model is different than that used in our single temperature 
model, i.e., fi - 2.0 instead of 1.5. However, this difference 
is expected, because, as already noticed in iDunne & EalesI 
(2001), single temperature models require lower values of 
than multi-component models. 

Constraints on y found in our study are in ve ry good 
agreement with estimates made bv lKovacs et al.l (l20IOh on local 
starbursts, i.e., y - 7.22 + 0.09. How e ver, u sing a sample of 
high-redshift starbursts, iKovacs et al.l (1201 Ol) found a lower 
value of y, i.e., y = 6.71 + 0.11. Thi s discrepancy migh t 
arise from the fact that to infer this value, IKovacs et al ] (120101) 
could only rely on uncertain MIPS-24 jim continuum estimates, 
extrapolated from broadband observations contaminated by 
PAH emission. 

We find an average emission diameter of 6 + 2 kpc 
(i.e., R = 3 kpc), which is consistent with estimates from 
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Fig. 3. Constraints on y, p and R obtained from a^- minimization analysis using 20 SMGs with PACS, SPIRE, submm and mm observations. 
These constraints correspond to our power-law temperature distribution model. Isocontours show the 99%, 95% and 68% confidence level. 



various studies using various high-resolution observations 
that have inferred diamet ers of order 1-10 kpc for SMGs 



(Chapman et al. 2004; Muxlow et al. 2005: 



Tacconi et al^'2006'. 



2008;.Bia2s & Ivison,200&; .Casev et aL.2009; lono et al. 2009, 


Lehnert et al. 


120091 CariUietal. 20101 Swinbank et al.1 poTq 


Tacconi et al. 


bold lYouneer et alJ 2010h. Kovacs et al. OOlol) 



found an emission diameter of ~2kpc for their high-redshift 
star-forming galaxies. As alr eady mentioned, t his discrepancy 
might arise from the fact that iKovacs et al.l (12010) relied on ex- 
trapolated MIPS-24 yum continuum measurements to make these 
estimates. We would like to stress that while our constraints on 
R are in line with previous estimates, its exact value should still 
be treated with caution. Indeed, robust constraints on the size of 
the emitting region would require the use of a complex radiative 
transfer model, taking into acc ount the geometry of the star - 
forming regions. For example, Chakrabarti & McKee' ("20081), 
using a self-consistent radiative transfer model and assuming a 
spherical geometry, found /?(;~10kpc. The agreement, within a 
factor 2-3, between our findings is encouraging in view of the 
approximations of our simple model. 

Based on these results, we conclude that fi, y and R 
can be considered as universal for these 19 SMGs. Nevertheless, 
how representative are these 19 SMGs of the full 61 SMG 
sample ? Using a KS analysis, we find that the redshift distri- 
bution of these two samples are fully compatible but that their 
infrared luminosity distributions are slightly different (only 
40% of chance of being drawn from the same distribution). The 
sample of 19 SMGs exhibits slightly higher infrared luminosi- 
ties than the full SMG sample, a median Lir of 6 x 10'^ Lq 



versus 4 x 10'^ Lq. These 19 SMGs are therefore not a perfect 
subsample of our full SMG sample. However, because these 
two samples are also far from being incompatible, we consider 
that the inferred values of fi, y and R are universal for our 61 
SMGs. This assumption is further supported by the fact that 
these parameters provide a good description of the far-infrared 
SED of the rest of our SMG sample (see Section l4.2.2b . 

4.2.2. Fitting tine full SMG sample 

We now fit the full SMG sample (including their PACS 70 yum 
and PACS 100//m detections) leaving and Mdust as the only 
free parameters of the model. Results of these fits are shown 
in Fig. lA.ll Table [13] summarizes the results inferred from 
these fits. Uncertainties are estimated using the distribution 
of Tc, Mdust and Ljr values that correspond to models with 
x''' < xlr^n ^- ^'^^ most of our SMGs this model provides 
(even with fixed fi, y and R parameters) a very good fit to our 
data points xi,^^~^ for Nqof~3). Almost all the highest 
values (> 25) correspond to the lensed-SMGs with relatively 
low infrared luminosities and high dust temperatures. This 
might suggest that for these galaxies fi, y and R are slightly 
different. Consequently, the infrared luminosities and dust 
temperatures inferred for these galaxies using our power-law 
temperature distribution model might be biased. These possible 
biases are discussed later on in this section. 

As for the single temperature model, we would like 
to verify that fits of SMGs with only PACS and submm obser- 
vations, or only SPIRE and submm observations, are not biased 
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Fig. 4. {Left) Dust temperatures inferred from the combination of PACS only (or SPIRE only) together with submm observations, compared with 
the reference values inferred using PACS, SPIRE and submm observations. These comparisons are for our power-law temperature distribution 
model. Symbols are the same as in Fig.|2] {Right) Same comparison but for the inferred infrared luminosities. The dust temperatures and infrared 
luminosities of galaxies can be reasonably inferred from their PACS -I- submm or their SPIRE+ submm observations alone using a temperature 
distribution model. 



compared to fits of SMGs with PACS, SPIRE and submm 
observations. Therefore, we compare the dust temperatures 
and infrared luminosities that one would infer using only the 
PACS (or SPIRE) and submm observations and our power-law 
temperature distribution model (with = 2.0, y - 7.3 and 
R - 3kpc), to that inferred using the combination of PACS, 
SPIRE and submm observations. This analysis is based on 
50 SMGs detected by both PACS and SPIRE and results are 
shown in Fig. |4] We find that the dust temperatures and infrared 
luminosities inferred using the combination of PACS (or SPIRE) 
and submm observations are in very good agreement with those 
inferred using the combination of PACS, SPIRE and submm ob- 
servation: o-[T^^l - Tl^^^] = 1 .2 K ( cr[rf <=f - r.sPiRE] = 2.3 K ) 
and o-[Lf//Lf^^^] ="o.lOdex ( cr[Lfrf/LSPif ^ 0.09 dex) 
This agreement is even better than that obtained in the case of 
our single temperature model. Consequently, estimates made on 
SMGs with only PACS or only SPIRE observations can be used 
with confidence. 

One of the main advantages of this power-law tem- 
perature distribution model is that it provides robust esti- 
mates of the total infrared luminosity (Lir[8 - lOOOjum]) 
of galaxies. The left panel of Fig. |5] shows the difference 
between the infrared luminosity extrapolated from a single 
temperature modified blackbody model and that inferred 
from our power-law temperature distribution model. We 
find a very tight correlation between those two estimates, 
Iog(Lj^"i''-T) ^ 0.84(+0.02) X log(Lj|;'^'""^) + 2.0(+0.2). 
However, we observe that the single dust temperature model 
systematically overestimates the luminosity of galaxies at 



high infrared luminosities and underestimates the luminosity 
of galaxies at low infrared luminosities. These discrepancies 
can be explained by the fact that in our single temperature 
model we were using a constant bolometric-correction term to 
convert LiR[40-120//m] into LiR[8-1000/vm], while its value 
changes with dust temperature (as with infrared luminosity, 
because there is a broad - Lir correlation; see the right 
panel of Fig. |9]l. For example, at high infrared luminosity 
(i.e., LiR > 3 X IO'^Lq) all galaxies have > 25 K. At 
these temperatures, the bolometric-coiTection term is, in our 
power-law temperature distribution model, of the order of ~ 1.5. 
The difference between our constant bolometric-correction term 
of 1.91 and this one, fully explains the observed discrepancies. 
This bias is illustrated by the green dotted-dashed line in the left 
panel of Fig.|5] 

The right panel of Fig.|5]shows the comparison between the 
dust temperature inferred using a single-temperature modified 
blackbody and that inferred using our power-law temperature 
distribution model. There is a tight correlation between these 
estimates and a very small dispersion. However, we can observe 
significant differences between these two estimates. indicates 
the temperature of the coldest dust component while Tdust mea- 
sures an effective dust temperature, therefore it is not surprising 
that Tdust yields values warmer than Tc. Some lensed-SMGs 
significantly deviate from this Tdust-T^c relation. These galaxies 
correspond to the ones with the largest x\^\ values, suggesting 
that in these systems (3, y and R might be slightly different. 
These dust temperatures are systematically shifted towards 
lower values while the corresponding infrared luminosities are 
not affected (see the left panel of Fig. |5]). Consequently, when 
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Fig. 5. {Left) Comparison of the infrared luminosities inferred using a power-law temperature distribution model with those inferred using a single 
dust temperature model. Symbols are the same as in Fig.|2] The black solid line shows the one-to-one relation. The green dotted-dashed line shows 
the bias introduced in our single dust temperature model by the use of a constant bolometric-correction term of 1.91 to convert Lir[40- 120 ;um] into 
iiR[8- 1000 yum]. To compute this line we measure Lir[40- 120yum] and Li[{[8- 1000/jm] on a power-law temperature template library normalized 
to reproduce the - Lir correlation (see the red dashed line in the right panel of Figure|9]l. We then plot on the j-axis 1.91 x Lir[40- 120 /jm] and 
on the y-axis Lir [8 - 1000 ;um]. (Right) Comparison of the dust temperatures inferred using a power-law temperature distribution model (T^) with 
those inferred using a single dust temperature model (rdust)- The red dashed line show a linear fit to the Tc-Ldust relation, = 0.6 x Tii,st + 3 K. 
Symbols are the same as in the left panel. Note that indicates the temperature of the coldest dust component of the multi-component model 
while Tdust measures an effective dust temperature. 



studying the - Lir plane, one has to keep in mind these slight 
shifts, or refer to the Tdust - ^ir plane which is not affected by 
this effect. 

In the rest of the paper we use the infrared luminosi- 
ties derived using the power-law temperature distribution 
model, unless stated otherwise. 



5. Results and Discussion 

5.1. The Infrared luminosity of SMGs 

The nature of SMGs has been greatly debated for more than a 
decade and in particular the reliability of their measured extreme 
SFRs. Indeed, while simulations of major mergers are able to 
reproduce such extreme SFRs, simulations in a cosmological 
context have had great dif ficulties accounting for the estimate d 
SFRs and number counts dBaueh et aI.ll2005UDave et al.ll2010h . 
Thus, the question remains: are the infrared luminosities of 
SMGs overestimated ? Thanks to Herschel observations we 
can now assess this question by measuring the true infrared 
luminosity of SMGs, studying their evolution as function of the 
redshift and testing the quality of pre-Herschel estimates based 
on monochromatic extrapolations. 



Figure |6] shows the infrared luminosities of SMGs as 
a function of their 850 fim flux densitie^ and their red- 
shifts. Our results unambiguously confirm the remarkably 
large infrared luminosities of SMGs. The vast majority ex- 
hibit infrared luminosity larger than 10'^ L©, and some even 



have Li 



IR 



> 



^Q. The first, second and third quartiles 
of our sample are lO'^^L©, 10'^'^ Lq and 10'^'^ Lq, respec- 
tively. These infrared luminosities correspond to SFRs of 
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for sources with no 850 fim observations we used extrapolations 



assuming f} 



2 i e S ^^'™P*^^'^^^'^ 



X (^,, 



„/S50)'* where 



lOOMoyr"', 4OOM0yr"' and 630MQyr"', respectively (using 
SFR [Mo yi--'] = 1 X 10"'°LiR [Lq], assuming a Chabrier 
IMF and no significant AGN contribution to the far-infrared 
luminosity). The existence of this large sample of star-forming 
galaxies with extreme infrared luminosities illustrates the strong 
evolution with redshift of the infrared galaxy population: in the 
local Universe such luminous infrared galaxies are very rare 
but their comoving spac e density increases by a factor ~400 
betwe en z ~ and z ~ 2 (iMagnelli et al.l 1201 it IChapman et al.l 
l2005h . Consequently, the characterization of the mechanisms 
triggering their starbursts becomes crucial in order to obtain a 
good census of the star formation history of the Universe. 

We observe a weak trend between S^so and Ljr (left 
panel of Fig. |6|. However, this coiTelation is likely driven by 
selection effects. Indeed, since submm observations at low lu- 
minosity are biased towards cold dust temperatures (see section 
13.81 ), they miss the bulk of the star-forming galaxy population 
at low and intermediate infrared luminosities. This missing 
population should have warm dust component s and therefore 



, is the (sub)mm wavelength at which the SMG has been detected. 



relatively faint 850 yum flux densities (see also Chapman et al 
2004; Casevetal. 2009; Magnelli et al. 2010; Chapman et ^ 
2010; Magdis et al. 2010). This hypothesis is strengthened by 
the position of some of the lensed SMGs, which give us a 
glimpse into the bulk of the population of galaxies with low 
infrared luminosities. The underlying 5 sso - ^ir relation cannot 
be probed using a submm-selected sample. 

Submm observations have the great advantage of being 
subject to negative A:-correction which makes an galaxy equally 
detectable in the submm over a very wide range of redshift. 
Therefore, one can expect the redshift distribution of submm 
galaxies to be relatively uniform if there were no strong evolu- 
tion of the underlying galaxy population. Instead, we observe a 
strong correlation between the infrared luminosities of galaxies 
and their redshifts (right panel of Fig. |6l). This trend can be 
explained by an evolution of the underlying galaxy population 
and by selection effects. The increase with redshift of the 
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S(850^m) [mJy] redshift 

Fig. 6. {Left) Infrared luminosities as function of the submm flux density. Blue squares represent SMGs situated in blank fields while green 
diamonds represent lensed-SMGs. OFRGs from lMagnelli et ^ ( 120 IQ) are presented with left red arrows. The solid and dashed lines show the 
linear fit to the Lir - 585o relation and the Icr envelope (Lir[Lo] = ]^oii.33±o.29 x Sgjg' [mJy]). Dotted lines show the Lir -Shq relation followed by 
single modified (fj = 1.5) blackbody functions at 20, 35 and 50 K. (Right) Infrared luminosities as function of the redshift. The symbols are same 
as in the left panel but OFRGs are represented by red filled circles. Blue dotted, red dashed and green dotted-dashed lines present the lower limit 
of the parameter space reachable using our deep radio (i.e., 20 /iJy), PACS 160 jim. (i.e., 3 mJy) and MIPS-24/im (i.e., 20yuJy) observations of the 
GOODS-N field, respectively. Note that in these figures galaxies with high^- value do not lie in a particular region of these plots but are rather 
randomly distributed. 



number of very luminous SMGs is due to the evolution of the 
infrared galaxy population and a volume effect: at high redshift, 
the comoving space^ density of luminous infrared galaxies is 
larger ( Magnelli et all i201ll: IChapman et alJ l2005h and the 
comoving volume probed by our survey increases. On the other 
hand, the lack of low luminosity galaxies at high redshift is 
quite surprising. Indeed, simply due to a volume effect, we 
would expect to see many more low luminosity galaxies at 
high redshift than at low redshift. We argue that this trend 
can be easily understood as a pure selection effect. Indeed, as 
illustrated in the right panel of Fig. |6] the depth of the deepest 
radio observations used to provide robust multi-wavelength 
counterpar ts creates th e low bound a ry in the Lir - z plane. 
iPope et al.1 (^006) and iBanerii et all (1201 ll) argue instead that 
this lack of low-luminosity galaxies at high redshift could be due 
to an evolution of their SEDs. To be missed by submm obser- 
vations, those galaxies should exhibit hotter dust temperatures 
than low redshift galaxies of the same luminosity. This seems to 
be incompatible with the modest evolution with redshift of the 
T'dust - ^iR re lation observ ed up to z ~ 2 (iHwang et al.ll2010l; 
IChapinet all EzOll: Marsde n'etaDl201 ih . 

Using our reference infrared luminosities (i.e., inferred 
from the power-law temperature distribution model) we can 
now test the quality of pie-Herschel estimates. One of the most 
common pie-He rschel monochromatic extrapolations was based 
on the MIPS-24 fim flux densities and the lCharv & Elbazl (1200 ll 
hereafter CEOl) SED library. We applied these extrapolations to 
our SMG sample and compared those estimates (hereafter Lj^) 
to our reference infrared luminosities (left panel of Fig.|7]i. 

Our results reveal that the use of the MIPS-24 fim emission 



and of the CEOl SED library yields inaccurate estimates 
of the infrared luminosities, characterized by a large scatter 
(cr[log(Lj^/L[^*^)]~0.47 dex) and a systematic overestimate for 
the most lumi nous galaxies. These results are in line with 
conclusions of 'Hainline et al.' (l 2009b s tt idying SMGs and of 
Pacovich et al. (2007), Murnh vet al] i200% . iNordon et all 



6201 a i2012l) and Elba z et al.l ' (l2010H201 l lfstudving bolomet- 
rically selected high-redshift galaxies. Our study also agrees 
with the fact that the overestimate of the infrared luminosity 
by the MIPS-24 ;um flux density and the CEOl SED library 
occurs at z > 1.5, i.e., when the MIPS-24 /im passband starts 
probing rest-frame wavelen gths dom inated by PAH emission 
(Nordon et al. 201^ |20T1 lEIbaz et al. 2010, l20Tlh . Indeed, 
SMGs with infrared luminosities below 10'^ are all at z < 1 .5 
and exhibit better agreement between L^^ and D^^. 

All these studies show that the SEDs of star-forming 
galaxies strongly evolve with redshift. This evolution might 
be interpreted as a modification of t he physical conditions 
prevailing in t heir s tar-forming regions. lElbaz et al.l (l20TTl) and 
iNordon et akl (12012 ) found that the SEDs of these high-redshift 
galaxies with extreme star-formation could be described using 
local SEDs o f less luminou s galaxies (see also iPapovich et al] 
.2007; ,MagnelIi et al.ll201 lb . This SED evolution is thus likely 
due to an increase of the PAH emission strength: the star-forming 
regions in those extreme high-redshift starbursts might be less 
compact than in their local a nalogues (i.e., ULIRGs), resulting 
in stronger PAH emission (IMenendez-Delmestre et al.l l2009l) . 
This hypothesis is supported by the observation s in SMGs of 
larger star-fo rming regions than in local ULIRGs dTacconi et al.l 
20061 120081 [2010) and by the observations of stron g PAH 
signatures in their IRS spectra (Lutz et al..,2005; .Valiante et al.l 
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Lir"™" [Lq] from CE01 + 24um Lir""™" [Lq] from 1.4GHz 

Fig. 7. Infrared luminosities for submm sources detected at 24 yum and 1.4 GHz. The x-axis shows the infrared luminosities extrapolated from the 
MIPS-24yum (left) or the radio (right) flux density, using the CEOl library or the FlR/radio correlation (with q = 2.34), respectively. The -y-axis 
shows the ratio of the infrared luminosities extrapolated from the MlPS-24 fim or radio flux density and the reference infrared luminosities inferred 
from our power-law temperature distribution model. The symbols are same as in Fig[6] 
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Fig. 8. Evolution of {q} as function of the infrared lumi nosity (left panel) and the redshift (right panel). On the left panel, solid and dashed lines 
show the local relation and its 1 cr dispersion as found b vlYun et alj i 200i|). On the right panel, solid and dashed lines show the redshift evolution 
of (q) oc (1 +^)-o.i5±o.o3 fjQjj^ local value as inferred inl lvison et al. Il l2010a l). The symbols are same as in Fig|6] Note that here, we did not attempt 
to correct for any incompleteness, e.g., using a Kaplan-Meier estimator, and biases introduced in our sample. So these results should be taken with 
caution because they only apply to our specific selection function, i.e., SMG with spectroscopic redshift estimates mainly obtained through robust 
radio identifications. 



120071: iPope etal]|2008l:lM"enendez -Delmestre et al ]|2007ll2009l) . olation was to use radio flux densities and the local FIR/radio 

correlation dHelou e t al. 1988; ,Yun et al.,.2001,) . 

Another popular pTe-Herschel monochromatic extrap- 



? = log 



dm 



3.75 X 1012 X Li.4GHz[WHz-i] 



(5) 



where LpiR is the infrared luminosity from rest frame 40 fim 
to 120 fim and Li.4ghz is the ^-corrected radio luminosity den- 
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Fig. 9. (Left) Dust temperature-luminosity relation inferred from our single temperature model. The symbols are same as in Fig [6] Red circles 
present the OFRG sa mple of [M agnelli et al. (2010). The striped area presents results for SMGs extrapolated by Cha pman et alJ (2005) from radio 
and submm data. The lCha pman et al. ( 2003) derivation of the median and interquartile range of the Tjust-AR relation observed at z ~ is shown 
by solid and dashed-do tted lines, linearly extrapolated to 10''' Lq. The dashed line represent the dust temperature-luminosity relation derived in 
iRoseboom et"^ ^20 11) for mm-selected sample observed with SPIRE and assuming a single modified blackbody model. {Right) Dust temperature- 
luminosity relation inferred from our power-law temperature distribution model. Symbols are the same as in the left panel. The red dashed line 
presents the T^-Li^^ relation inferred from a least-square second degree polynomial fit. 



sity (h ere we assume a standard radio slope a - 0.8; llbar et all 
120101) . In the following, we derived the infrared luminosities of 
our galaxies using this FIR/radio correlation and {q ) - 2.34, as 
observed in the local Universe bv lYun etalJ(l2001h . Those esti- 
mates are compared to our reference values in the right panel of 
Fig.E] 

We find a tighter correlation between our reference infrared 
luminosities and those inferred using radio flux densities and the 
local FIR/radio correlation (cr[log(Lfjj"*7L5^f)]~0.29 dex). The 
accuracy of these extrapolations is also supported by the good 
agreement found between those estimates in our lensed SMG 
sample. Nevertheless, we also observe a trend with the infrared 
luminosity: at high luminosities, the FIR/radio correlation sys- 
tematically overestimates the luminosity. Since there is a tight 
correlation between Lir and z, one can suspect this trend to be 
driven by an evolution of {q) with redshift. As illustrated in the 
right panel of Fig. [8] this trend is indeed in very good agree- 
ment wi th the evolut i on of jq) proportional to (1 + ^y^'' '^^^^-^^ 
found in llvison et a n (l2010ai) . Nevertheless, one has to keep in 
mind that our sample cannot be used to probe the evolution of 
(q) with redshift, since it is, by construction via the radio identi- 
fications, biased towards galaxies with high radio flux densities. 
Therefore, because here we did not attempt to correct for any of 
these incompleteness, e.g., using a Kaplan-Meier estimator, any 
of our results on (q) should be taken with caution. The evolution 
of (q) could only been studied through carefully selected sam- 
ples and using radio stacking. So far, no clear conclusion on the 
evolution of ( q) with redshift has been made (see Ivison et al. 
I2010all20105 Roseboom et al. 1201 ID . 

5.2. The Tdust - ^ir plane 

The left panel of Fig. |9] shows the Tdust - i'lR plane inferred 
from our single temperature model. The use of this simple 
model provides a comparison to other studies. Compared 
to previous Hersch el-hased results (Magnelli et al. 2010; 
IChapman et aI.ll20Toh . our large SMG sample populates the low 
(i.e., LiR < IO'^ ^Lq) and high (i.e., 10'^ Lq > Lir) luminosity 



regions of the Tdust - ^ir diagram. This large dynamic range 
allows a clear characterization of the Tdust - i^iR correlation. 

The left panel of Fig. |9] clearly confirms the selection bias 
introduced by submm observations: At low luminosities SMGs 
are biased towards cold dust temperatures. The upper envelope 
of the SMG Tdust - ^ir distribution only depends on the depth 
of the submm observations (see Section ITSl ). The existence of 
a population of dusty star-forming galaxies missed by submm 
observations is corroborated by the presence, in the upper part 
of the Tdust - LiR diagram, of some of the lensed SMGs and the 
optically faint radio galaxies (OFRGs, iMagnelli et alj I20T0I ; 
Casey et al. 2009). 

Our SMG sample, tog ether with our lense d SMG sample 
and the OFRG sample of .Magnelli et alj (I201Q) . suggests that 
high-redshift dusty star-formin p galaxies e xhibit a wide range 
of dust temperatures (see also 'Casev et al.l 120091 : iMagdis et al.l 
2010). This might indicate that the Tdust - ^ir relation at high 
redshift has a higher scatter than locally. However, this conclu- 
sion can be driven by selection effects, because a significant 
fraction of the galaxies with intermediate dust properties are 
missed by our current sample. This missing populati on will 
probab ly reconcile our finding with those of Hwang Tt al.l 
(I2010I) . who found modest changes in the Tdust-iiR relation as 
function of the redshift using an LjR-selected sample of galaxies 
observed with Herschel. This conclusion is also strengthened by 
the fact that at high luminosities (i.e., few times 10'^ L©, where 
SMGs are a representative sample of the entire high luminosity 
galaxy population) SMGs exhibit dust temperatures that are 
in line with the Tdust - ^ir relation extrapol ated from local 
observations of Chapman et al. (2003) (see also IClements et aP 
I20 1 Ol: |Planck"Colla boration et al.ll2011h . 

As illustrated by the striped region in the left panel of Fig. 
|9] our dust temperatures and infrared luminosities largely agree 
with those extrapolated by pie-Herschel studies using the local 
FIR/radio correlation. This agreement of course reflects the 
broad consistency found between the local value of (q) and that 
observed in our sam ple (see Sectionl5.ll). Ou r results also agree 
with those found by iRoseboom et al.l (1201 Ih on a mm-selected 
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sample observed with SPIRE and assuming a single modified 
blackbody model (see the dashed line in the left panel of Fig.|9]l- 

From the wide range of dust temperatures, we can 
conclude that although submm observations are very useful to 
select extreme star-forming galaxies at high redshift, they cannot 
be used to obtain a complete census of the dusty star-forming 
galaxy population with relatively low infrared luminosities 
(LiR < 10'^'^ Lq). This census is now possible using bolomet- 
ric selections prov ided by deep Herschel observations (e.g., 
iMagdis et al. l l2mol) but still limited to relatively low redshift 
galaxies (z < 2.5) due to the positive ^-correction affecting 
Herschel data. In the near future, very deep mm observations 
provided by the Atacama Large Millimeter Array (ALMA) 
might help to obtain this census even at high redshift. 

The right panel of Fig. |9] shows the Tc - Lir plane in- 
ferred from our power-law temperature distribution model. This 
plane cannot be compared to any pie-Herschel studies. is the 
temperature of the coldest dust component while Tdust gives an 
average dust temperature. Thus, is systematically lower than 
T^dust, but their relative variations are tightly correlated (see also 
Fig. |5]l. Conclusions that one can draw from the - Lir plane 
are the same as those drawn from the Tdust - ^ir plane. 

We fitted the - log(LiR) and Tdust - log(iiR) relation 
with a second order polynomial function and studied the scatter 
around these fits. We find <tt^ = 1.9K and ctt-^^ ^ - 3.8 K. The 
decrease of the scatter (by a factor 2) is in line with expectations 
from the relation between and Tdust, i-e., a factor 1.7 because 
Jc = 0.6 X Tdust + 3 K. We note that our single-temperature 
model is also sensitive to the rest-frame wavelengths used in 
the fits; even if all galaxies at a given infrared luminosity have 
had the same dust temperature, our single-temperature model 
would still be affected by their redshift distribution, i.e., by the 
rest-frame wavelength probed by the PACS 160/im data point. 
This redshift distribution would then introduce an artificial 
Tdust scatter. In contrast, our power-law temperature distribution 
model is less affected by this effect, since it is constructed to 
reproduce cold and warm dust components. 

5.3. The spectral energy distribution of SMGs 



As discussed in section 15.21 SMGs can not be treated as a 
homogenous galaxy population, because they probe wide ranges 
in infrared luminosity and dust temperature. Moreover, while 
at high infrared luminosity SMGs are a representative sample 
of the underlying high luminosity galaxy population, at lower 
infrared luminosities, SMGs are only a subsample of the entire 
infrared galaxy population, and are biased towards cold dust 
temperatures. Thus, the SEDs of SMGs have to be analysed 
as a function of their infrared luminosities. Figure [TOl presents 
the photometry of our SMGs split into four different infrared 
luminosity bins, i.e., Lir < 10'- Lq, 10'^ Lq < Lir < 10'^-^ Lq, 
10'^ ^Lq < Lir < 10'^ Lq and 10'^ Lq < Lir. In these panels, 
we show the mean SED inferred from our power-law temper- 
ature distribution model. These SEDs correspond to p - 2.0, 
y - 7.3, R = 3 kpc and to the mean Mdust and inferred for 
the galaxies of the bin. In these panels, we also show the CEOl 
SED corresponding to the mean infrared luminosity inferred for 
the galaxies of the bin, i.e., the CEOl SEDs are not fitted to the 
photometry of the individual galaxies here. 

At high infrared luminosities, the peak of the CEOl 
SED template is in agreement with that of our power-law 







J 10'° 








Teslframe 

Fig. 10. Mean rest-frame SED of SMGs for four infrared luminosity 
bins, from bottom to top: Lir < 10'^ Lq; 10'- Lq < Lir < lO'-^ Lq; 
10'^ ' Lq < L,R < 10'-'' Lo and 10" Lq < Lir. The solid lines show 
the power-law temperature distribution SED corresponding to the mean 
dust mass and dust temperature of the bin. The photometry of each of 
the sources was slightly renormalized to match these SED templates at 
submm wavelengths. Dashed lines represent the CEOl template corre- 
sponding to the mean infrared luminosity of the bin, i.e., these templates 
were not fitted to the photometry of our SMGs. 



temperature distribution model. This indicates that in this range 
of luminosities, the local Ldust - Lir relation used in the CEOl 
library does not significantly evolve with redshift. In contrast, 
the MlPS-24yum observations are systematically above pre- 
dictions from the CEOl SED template. As already mentioned, 
these discrepancies are likely due to an increase of the PAH 
emission strength in these galaxies and produce inaccurate 
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infrared luminosity extrapolati ons from the MIPS-24 //m flux 
densi ty using the CEOl library (lElbaz et al.ll201 it iNordon et al] 
l20T2h . 

As we go to lower infrared luminosities we observe larger 
discrepancies between the peak of the CEOl SED and that of 
our power-law temperature distribution model. While the CEOl 
SED templates follow the local rdust - ^ir relation, our SMG 
sample is more and more biased towards cold dust temperatures. 
At such low infrared luminosities, the SMG population repre- 
sents the low-temperature-end of the real Tdust - i'lR distribution 
(see Fig. |9l). In this low luminosity range, we note the better 
agreement than at high luminosities between observed and 
predicted MIPS-24 fim. 



6. Toward a better understanding of the nature of 
SMGs? 

Our results unambiguously reveal the diversity of the SMG 
population. Some of these galaxies exhibit extreme infrared 
luminosities, with no local analogues (Lir > lO'-'L©), 
while others have relatively low infrared luminosities 
(10^^ Lq < Lir < 10'^ Lq). Is this diversity reflecting dif- 
ferences in the mechanisms triggering their SFRs? 

Recent hydrodynamic simulations, coupled with ra- 
diative transfer calculations, have found that while SMGs 
with relatively low infrared luminosities can be created by 
different scenarios (two gas rich galaxies soon to merge and 
observed as one submm source, or an isolated star-forming 
galaxy with large gas fraction), SMGs with the most extreme 
infrared luminosities/SFRs (i.e., Lir > 10'^^ Lq, equivalently 
SOOMoyr"') can only be indu ced by strong s t arburs ts at the 



coalescence of major mergers dHavward et alj 201 Ih . These 
results are consistent with those of Dave et alT i 20 lo') who found 
that SFRs induced by a secular mode of star formation reach at 
most, at z ~ 2, a value of -SOOMoyr"^ (i.e., Lir ~ 10'^ '' Lq). 
This value of ~500Moyr ' can thus be considered as the 
maximum non-merger SFR" (hereafter SFR;^™'^^ and be 
used to separate, at z ~ 2, merger-induced starbursts from 
galaxies with a secular mode of star formation. Moreover, in a 
steady-state between SFR and gas accretion, one could expect 
SFR^^"'"' and the ga s fraction of galaxi es to be strongly related 
(iBouche et alJ 120101: iDave et al.ll20Tlh . Therefore, SFR^™'^' 
should decrease at low redshift with the gas fraction of galaxies 
(Tacconi et al. 2010; Geachetal. 2011). Qualitatively this 
assumption is supported by observations of local ULIRGs, 
which are mostly associated with major mergers but which 
exhibit SFRs lower than ~ 500 M© yr"\ likely because they have 
relativ ely low gas fraction (~ 10%; see Fig. 9 of Saintonge et al. 
Hon]). Therefore, in the redshift range z - 0-2, we can separate 
merger-induced starbursts from non major-merging ones using 
a threshold of 500 x (1 -H z)?^^ yi""', while at z > 2, we can 
use a threshold of 500Moyr Here, the redshift dependence, 
(1 + 2)2^2 = ((1 + z)/3)^'^, cornes from the evolution of the gas 
fraction found in iGeach et"ai] (1201 lb . Using a Chabrier IMF, 
these SFR thresholds correspond to the most luminous SMGs 



of our sample, i.e., Lir > 10^^ '' Lq x (1 -H z)^l"; at < z < 2 and 



2.2 



z=2 



> 10'^ ' L0 atz > 2. 
A correlation between the SFR and the stellar mass of 
star-forming galaxies has been observed over the last lOGyr 
of lookback time (SF Roc or SSFR^SFR/M,oc Mr' 

wifli 0.5 < a < l -O; iNoeske et all 1200 7; ElbazetalJ 



iDaddi et all l2007bt iRodighiero et all 120101: lOhver et al 



2007 



201c ; 



iKarim et all 1201 it iMancini etani201lh . The existence of this 
"main sequence of star formation" (MS) is usually interpreted 
as a piece of evidence that the bulk of the star-forming galaxy 
population is forming stars gradually with long duty cycles. 
Galaxies situated on the main sequence would be consistent 
with a secular mode of star formation, likely sustained by a 
continuous gas accretion from the IGM and along the cosmic 
web (Dekeletal. 2009; Dave et al. 2010), while star-forming 
galaxies located far above the main sequence would be con- 
sistent with strong starbursts with short duty-cycles, mainly 
triggered by major mergers. In that picture, to separate galaxies 
triggered by major mergers from those with secular mode of 
star formation, one should use the offset of a galaxy with respect 
to the MS, r ather th an simply u sing i ts infrared lurn i nosity 
(IWuvts et al.' I20T1I: lElbaz et all I201TI iNordon et alj lIoTa 
[Rodighiero et al.ll201 ll Magnelli et al., in prep.). 

There are thus two ways to identify major-merger 
induced starbursts. In the following, we apply these two 
criteria to our SMG sample, compare their results, and 
more importantly test their ability to effectively select 
major-merger induced starbursts. For the criterion using 
the offset of a galaxy with respect to the main sequence (i.e., 
Alog(SSFR)MS = log[SSFR(galaxy)/SSFRMs(M.,z)]), we use 
the stellar masses of 39 blank field SM Gs derived in Sec t tonl3.7l 
and the definition of the MS given by Rodighi ero et al ] (1201 Oh . 
i.e., log(SSFR)MS = Q'log(M.)+je where (a , ^S) = (-0.27, 2.6), 
(-0.51,5.3) and (-0.49,5.2) at 0.5 <z < 1.0, 1.0 <z< 1.5 and 
z > 1.5 , respectively. We adopt the definition of Rodighiero et all 
(l20Tol) for consistent use of the FIR as a star-formation indicator. 
None of our results strongly depend on this specific definition. 

In the left panel of Fig. [TT] we observe that our SMGs 
are systematically above the main sequen ce of star-formation , 
consistently with p r evious findings (e.g. , iDaddi et aTl l2007bt 
iHainUne et al.ll201 ll; l^dlow et al.l l201 iV Nevertheless, while 
low luminosity SMGs are within Icr from the MS, SMGs above 
our merger-induced starburst separation (i.e., SFR^^™'^') are at 
least 2cr above it. This segregation shows that for the relatively 
narrow range of stellar masses probed by our SMG sample 
(1 X 10'° M. - 5 X 10" M.), a simple cut in SFR allows us to 
accurately select the galaxies lying above the main sequence. 
The fact that these two independent criteria (one is based 
on hydrodynamic simulations while the other is empirically 
derived using duty cycle arguments) select the same sample of 
galaxies strengthen their accuracy and therefore supports the 
assumption of a major-merger induced scenario. We note that 
the SFR/luminosity criterion selects galaxies located ~ldex 
above the main sequence. This is consistent with values used 
in studies selecting merger-induced starbur sts based on their 
location with respect to the main sequence (lElbaz etal J I20TTI: 
INordon et all 1201 21: iRodighiero et allboilh . We conclude that 
for our specific SMGs sample these two criteria are equivalent. 

Half of the galaxies in our sample (29 SMGs) have 
SFRs above our merger-induced starburst separation (i.e., 
SFR^™'^'; hereafter we call these galaxies luminous-SMGs, 
because they have Lir > IO'^-^Lq X (1 -h zfj2 at < z < 2 
and Lir > 10'^'^ Lq at z > 2). Their median infrared luminosity 
is 6.4 X 10'^ L0, their median is 27 K and they are at least 
2cr above the MS of star formation. The high dust temperatures 
of these luminous-SMGs agree with those observed in local 
ULIRGs (see the agreement between the mean SED of these 
SMGs and the CEOl template in the top panel of Fig.fTO]). This 
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Fig. 11. (Left) Distribution of "distance" witii respect to tiie main sequence observed in our SMG sample having accurate stellar masses estimates 
(empty histogram). The hatched histogram shows the distribution observed in a subsample of luminous SMGs with SFR>SFRj^';!;"'"', i.e., Lir > 
10'^-^ Lo X (1 + z)l'^2 at < z < 2 and Lir > lO'"-^ Lq at z > 2. (Right) Dust temperature of SMGs as function of their distance with respect to the 
main sequence of star-formation. Blue squares show luminous SMGs with Lir > 10'^ ' Lq x (1 + z)^^2 at < z < 2 and Lir > 10'--^ Lg at z > 2. 
Green squares show SMGs with infrared luminosities below these thresholds. Red points represent the OFRGs, i.e., galaxies with relativ ely low 
infrare d luminosities below our treshold. In both plots the location of the main sequence as function of the redshift is taken from .Rodighiero et al.l 
( 120101) . The 1 cr scatter around this main sequence is illustrated by the shaded area. 



agreement suggests similar physical conditions prevailing in the 
star-forming regions of local ULIRGs and those of luminous- 
SMGs. Since local ULIRGs are triggered by major mergers, this 
suggest that luminous-SMGs might also be produced by major 
mergers. 

The relatively high dust temperatures of the luminous-SMG 
subsample (compared to the rest of the SMG population, see the 
right panel of Fig. fTTT l also agrees, qualitatively, with the large 
increa se of the dust temperature predicted by iHavward et al.l 
(1201 Ih at the coalescence of their major merger simulations. 
To quantitatively confirm this agreement we compare our dust 
temperatures with th ose of the hydrodynamic simulations of 
iHavward et al.l(l201 ih . First, we redshifted their simulated SEDs 
to match our SMG redshift distribution, second, we convolved 
these SEDs with the PACS, SPIRE, submm and mm filters and, 
third, we applied cuts in flux densities to match the properties 
of the GOODS-N field (i.e., a field with deep submm and 
Herschel observations, probing a large dynamic range in the 
Tc - iiR plane). Then, we fitted our power-law temperature 
distribution model with y6 = 2.0, y - 13 and R - 3> kpc to 
this set of simulated SEDs, leaving Mjust and as the only 
free parameters of the modeCB As for our data, the power-law 
temperature distribution model provides a good fit to the 
simulated SEDs, characterised by reasonably low values (i.e.. 



" If we constrained fi, y and R on the simulated SEDs, we find 
j8 = 1.6 ± 0.2, 7 = 8.7 ± 0.7 sad R = 2 ± Ikpc. These values are 
different that those obtained on our SMGs and lead, systematically, to 
higher dust temperatures (AL^ ~ 7 K). Nevertheless, we believe that 
using these constraints will not provide a fair dust temperature com- 
parison with our SMGs. First, while the exact values of f}, y and R 
strongly affect the inferred T^, the location of the FIR peak of the simu- 
lated SEDs stays unchanged. Simulated SEDs of major mergers peak at 
shorter wavelength than those of isolated starburst, and the localization 
of these peaks are consistent with those of our SMGs. Second, if the 
constraints on y and R from our SMGs do not provide the optimal fit 
to the simulated SEDs, they still provide a fairly good fit to them. Third, 
the simulated SEDs cannot be used to constrain yS, y and R because they 
do not yet include stochastically heated very small grains (Hayward et 
al., in prep). 



~ 8 for Ndof ~ 3). We find that simulated galaxies populate the 
same region of the Tc - Lir plane as our SMG sample. Extreme 
infrared luminosities (i.e., Lir > 10'^^ L©) are indeed only 
observed in simulat ions of strong starbur sts at the coalescence 
of major mergers (IHavward etalJ 1201 Ih . Simulations of two 
gas rich galaxies soon to merge (i.e., at an epoch where tidal 
effects have not yet caused strong starbursts) and observed as 
one submm source, always have lower infrared luminosities 
(i.e., Lir < 10'^'^ L©). In the simulations, strong starbursts at the 
coalescence of major mergers exhibit higher dust temperatures 
(Jc ~ 28 K) than isolated starbursts (T^ ~ 22 K). The agreement 
between the dust temperatures of simulated major mergers 
and that of our luminous-SMGs (i.e., is 27 K) supports 
the assumption that these luminous-SMGs are observed at a 
late-stage of a major merger. 

We conclude that the most luminous SMGs exhibit 
properties, including their extreme infrared luminosity 
(Lir = 6.4 X 10'- Lq), their hot dust temperature (rc = 27K) 
and their location with respect to the main sequence (> 2cr), 
which favour the scenario in which they correspond to intense 
starbursts with short duty-cycles, mainly triggered by major 
mergers. On the other hand, SMGs with low infrared lumi- 
nosities exhibit properties, including their relatively cold dust 
temperatures (T^ - 20 K) and their location with respect to 
the main sequence (within ~2cr), which favour the scenario 
of isolated star-forming galaxies or pairs about to merge, i.e., 
at a time where tidal effects have not yet caused strong star- 
bursts. The distinction between these two modes (i.e., isolated 
star-forming galaxy or early-stage major merger) cannot be 
assessed using our data. However, the existence in some SMGs 
of two galaxies about to merge and being contained in the 
same submm bea m is confirmed by so me high resolution CO 
line observations (iTacconi et al.1 l2008l HDF242 aka GN19) 
and submm continuum obser vatio ns (lYounger et al.l 12009 : 
Kovicsetalj[2010; Wang et al. 200). 

We stress that the conclusions drawn for the low luminosity 
SMG population should not be extrapolated to the entire low 
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luminosity galaxy population. Indeed, at these luminosities 
SMGs do not constitute a representative sample of the under- 
lying population (see Section [378T l. In particular, galaxies with 
relatively low infrared luminosities but warm du st (e.g., the 
OFRG s) might be triggered by major mergers dCasev et alj 
I20T1I) . 

In the right panel of Fig. [TT] we observe a clear corre- 
lation between the location of a galaxy with respect to the 
main sequence (i.e., Alog(SSFR)Ms) and its dust temperature 
Tc. Nevertheless, because our sample is affected by strong 
selection biases (in term of luminosity and dust temperature 
as well as in term of being preferentially optically bright), this 
Tc - Alog(SSFR)MS relation has to be treated with caution. 
At low infrared luminosities our sample is biased towards 
galaxies with cooler dust temperatures. Thus, we can expect 
"main sequence" galaxies to exhibit a broader range of dust 
temperatures than our current SMGs sample, i.e., weakening 
the Tc - Alog(SSFR)MS relation. On the other hand, the location 
of the OFRGs in this figure seems to qualitatively confirm the 
existence of a T^- Alog(SSFR)MS relation. Indeed, even at 
relatively low infrared luminosities (10^^ < Lir <, IO'-^Lq) 
and over the same range of redshift (i.e., 1 .0 < z < 2.5), galaxies 
with hotter dust temperatures (i.e., the OFRGs) are more offset 
from the MS than galaxies of the same infrared luminosities but 
with cooler dust temperature (i.e., the SMGs represented with 
green squares in the right panel of Fig. fTTT i. The existence of a 
Tc - Alog(S SFR)MS rel ation is also consisten t with the fact that 
lElbaz et al.l (|2011) and iNordon et al.1 (1201 2h find a coiTelation 
between Alog(SSFR)MS and the SED properties of star-forming 
galaxies. The PAH-to-LjR ratio of main sequence galaxies is 
constant, but decr eases with i n creasi ng of fset above th e main 
sequence. Finally, lElbaz et all ( 1201 lb and IWuvts et all (1201 Ih 
also find that Alog(SSFR)MS correlates with the compactness 
of the star-forming region; the SFR density of main sequence 
galaxies is roughly constant while it increases with increasing 
offset above the main sequence. All these correlations are strong 
observational support of the physical interpretation given to the 
main sequence of star-formation. Galaxies offset from the main 
sequence, likely triggered by major mergers, have compact 
star-forming regions resulting in warmer dust temperatures and 
weaker PAH emission. 

7. Summary 

Using the Herschel PACS and SPIRE observations of several 
deep cosmologcial fields, we study in detail the far-infrared 
properties of a sample of 61 SMGs which have secure redshift 
estimates. We find that at high infrared luminosities this sam- 
ple provides a good representation of the entire SMG population 
and more generally of the entire high luminosities star-forming 
galaxy population. At low infrared luminosities, our sample is 
less representative, because it is biased towards low redshift 
galaxies with cooler dust. Dust properties of these SMGs are in- 
ferred using two different approaches. First, we use a single dust 
temperature modified blackbody model which provides a very 
simple description of the dust emission of galaxies and allows 
comparisons with all pie-Herscliel estimates. Then, in order to 
obtain a better description of the Wien side of the dust emission, 
we use a power-law temperature distribution model. This model 
provides an accurate description of the rest frame far-infrared 
SEDs of SMGs. From this model we can constrain the dust emis- 
sivity spectral index, the characteristic emission diameter, the 
temperature index, the dust temperatures and the infrared lumi- 



nosities of SMGs. These properties are analysed and put into 
perspective with the more general question of the formation and 
evolution of star-forming galaxies. Our main conclusions are: 

1 . We find that a single dust temperature model provides a good 
description of the far-infrared peak and Rayleigh-Jeans part 
of SED of SMGs, but fails to reproduce its Wien-side. The 
dust temperatures and infrared luminosities inferred using 
the combination of only PACS (or only SPIRE) and submm 
observations are in very good agreement with the reference 
estimates based on PACS-(-SPIRE+submm data. 

2. Using a power-law temperature distribution model we ob- 
tain a good description of the far-infrared SED of SMGs at 
its peak, on the Rayleigh-Jeans side and on the Wien side. 
Using this model and the combination of PACS, SPIRE and 
submm observations, we obtain constraints on the dust emis- 
sivity spectral index of SMGs y6 = 2.0 + 0.2 and the temper- 
ature index y = 7.3 + 0.3. The dust emissivity spectral index 
found in our sample is in line with estimates by Dunne & 
Eales (2001). 

3. We find that luminosity extrapolations based on the radio 
emission are considerably more reliable than those based on 
the mid-infrared emission and the Charv & Elbaz (2001) li- 
brary. For our sample, the FIR/radio correlation is parame- 
terized with (q) - 2.0 + 0.3. However, this value could not 
be applied to the full high-redshift star-forming galaxy pop- 
ulations because our sample is not well-suited to study the 
evolution of (q) with redshift. 

4. Our study unambiguously reveals the diversity of the SMG 
population, which probes large ranges in infrared luminosity 
(from Lir ~ 2 X 10" Lq to ~ 3 x 10'^ Lq) and dust temper- 
ature (from Tc = 14 K to - 36 K) and is strongly bi- 
ased towards galaxies with cold dust. This bias decreases at 
high luminosities, and at Lir > lO'^^L©, SMGs are a repre- 
sentative sample of the entire high infrared luminosity star- 
forming galaxies population. At lower infrared luminosities, 
a complete census on the high-redshift star-forming galaxy 
population requires the use of the bolometric selection pro- 
vided by deep Herschel observations. 

5. Our study clearly reveals that some SMGs exhibit extreme 
infrared luminosities (Lir > lO'^^ L©) which correspond to 
SFRs of > 500 Moyr"' . We also observe that these luminous- 
SMGs exhibit warm dust temperatures (rc=27K) and are 
outliers of the main sequence of star-formation (~2cr above 
it). The extreme SFRs of these luminous-SMGs are diffi- 
cult to reconcile with a secular mode of star formation (e.g., 
iDave et ani2010h and could correspond to a merger-driven 
stage in the evolution of these galaxies. This hypothesis is 
supported by the fact that these SMGs exhibit warm dust 
temperatures consistent with estimates from hydrodynamic 
simulations of maj or mergers coupl ed with radiative trans- 
fer calculation (Ha vward et al.ll201lh . and that as outliers of 
the main sequence they are commonly assumed to be intense 
starbursts with short duty-cycles, likely triggered by major 
mergers. 

6. At low infrared luminosities, the dust temperatures and the 
infrared luminosities of SMGs are consistent with a secular 
mode of star formation. This hypothesis is also supported 
by the fact that those galaxies are situated close the main 
sequence of star-formation and hence are assumed to have 
large duty-cycles of star formation. 
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Table 1. Main properties of the PEP/HerMES fields used in this study. 





Field 




70 yum 


100 iim 




160 yum 




250 


350 yum 




500 yum 


Eff. Area 3o- 
arcmin- mJy 


Eff. Area 
arcmin^ 


30- 
mjy 


Eff. Area 
arcmin- 


3^" 
mJy 


EfF. Area 
arc mi 


mJy 


Eff. Area 
arcmin^ 


30-" 
mJy 


Eff. Area 30-" 
arcmin- mJy 


GOODS-S 


03''32'", 


-27°48' 


200 1.1 


200 


1.2 


200 


2.4 


400 


7.8 


400 


9.5 


400 


12.1 


GOODS-N 


12''36'", 


+62° 14' 




200 


3.0 


200 


5.7 


900 


9.2 


900 


12 


900 


12.1 


LH 


10''52'", 


+57°28' 




648 


3.6 


648 


7.5 


900 


11.5 


900 


16.8 


900 


24.3 


COSMOS 


00''00"', 


+02° 12' 




7344 


5.0 


7344 


10.2 


7225 


8.1 


7225 


10.7 


7225 


15.4 


A2218 


16''35"\ 


+66° 12' 




16 


2.6 


16 


6.2 


380 


18.0 


380 


9.0 


380 


9.9 


A1835 


14'"01'", 


+02°52' 




16 


3.4 


16 


6.9 


280 


18.0 


280 


11.5 


280 


12.2 


A2219 


leHo'", 


+46°42' 




16 


3.1 


16 


7.2 


280 


18.0 


280 


10.9 


280 


19.7 


MS 1054 


10''57'^ 


-03°37' 




50 


4.2 


50 


8.6 


500 


19.5 


500 


13.6 


500 


14.2 


CL0024 


00''26'", 


+ 17°09' 




36 


3.0 


36 


6.4 


280 


21.6 


280 


13.6 


280 


14.3 


MS0451 


04'" 54'", 


-03°01' 




16 


3.2 


16 


6.0 


280 


22.6 


280 


13.2 


280 


16.7 


A2390 


21''53'", 


+ 17°41' 




16 


4.0 


16 


7.9 


280 


38.8 


280 


20.0 


280 


20.3 


A370 


02''39'", 


-01°34' 




16 


3.0 


16 


6.0 


280 


18.0 


280 


12.0 


280 


13.1 


A1689 


13'" 11'", 


-01° 20' 




16 


1.8 


16 


4.0 


280 


21.0 


280 


12.4 


280 


13.2 



Notes. 

" in deep 160, 250, 350 and 500 /im observations, r.m.s values include confusion noise. 



Table 2. Main properties of the ancillary used in this study. 



Field 


24 


yum 


SCUBA-850yum 


LABOCA-870yum 


AzTEC-l.lmm 


MAMB0-1.2mm 


Radio 




3a- 


Ref." 


30- 


Ref." 


3o- 


Ref." 


3cr 


Ref." 


3o- 


Ref." 


30r 


Ref." 




fj]y 




mjy 




mjy 




mJy 




mjy 




^Jy 




GOODS-S 


20 


(1) 




N/A ... 


3.6 


(2) 


1.8 


(3) 




. N/A . . . 


20 


(4) 


GOODS-N 


20 


(1) 


3-12 


(5) 




N/A ... 


3 


(6) 


3 


(7) 


15 


(8) 


LH 


30 


(9) 


6 


(10) 




N/A ... 


2-4 


(11) 


2-3 


(12) 


16 


(13) 


COSMOS 


45 


(14) 




N/A ... 


16-30 


(15) 


3.9 


(16) 


3 


(17) 


34 


(18) 


A2218 


240 


(19) 


4.5 


(20) 




N/A ... 




N/A ... 




N/A ... 




N/A ... 


A1835 


50 


(19) 


5.1 


(21) 




N/A ... 




N/A ... 




N/A ... 


48 


(21) 


A2219 


30 


(19) 


5.1 


(22) 




N/A ... 




N/A ... 




N/A ... 


120 


(22) 


MS 1054 


36 


(19) 


4.5 


(23) 




N/A ... 




N/A ... 




N/A ... 


200 


(24) 


CL0024 


54 


(19) 


4.5 


(25) 




N/A ... 




N/A ... 




N/A ... 


45 


(25) 


MS0451 


42 


(19) 


12 


(22) 




N/A ... 




N/A ... 




N/A ... 


200 


(24) 


A2390 


30 


(19) 


6.6 


(25) 




N/A ... 




N/A ... 




N/A ... 


300 


(25) 


A370 


30 


(19) 


5.7 


(25) 




N/A ... 




N/A ... 




N/A ... 


30 


(26) 


A1689 


42 


(19) 


3.0 


(23) 




N/A ... 




N/A ... 




N/A ... 




N/A ... 



Notes. 

references: 



(7) fereve et al." (2008'), (8) 



(1) Magnelli etaD fcOllh . (2) IWeil3 et al. (2009b), (3) Scott et al. (20iU), (4) 'Miller et al." ('2008^, (5) 'Borvs et al. (2003|), (6) [P erera et al. (2008 

iMorrison etal] t2010.) , (9) Egami et al. (in prep.), ( 10).Coppin et al. (2006), (11) Austermann et al. (2010), (12) .Greve et al. (2004). (13) lBiggs & Ivison, (i200l„ , . . 

(15) Albrecht et al. (in prep .), (1 6)IScott et alj (|2008|), (17) Bertoldi et al. (2007), (18) Schinnerer et al. (20l3), (19 ) Valchanov et al. (i n prep.), (20) lKneib et al.l ( l2004h . (21) llvison et'al] OOOO 
(22) IChapman et alj ( 120021) . (23) lKnudsen et alj ( 120081) . (24) iStocke et al.Ul999l) . (25) iSmail et al.j 82002) and (26) IIvison et"!!] ( 119981) . 



I fere 

_(14) lLe Floc'h et al.. C2009 



•9 



o 



Table 3. Mid- and far-infrared properties of our GOODS-N SMG sample. 



Reference submm source and its counterpart PEP/HerMES muW-wavelength counterpart 



name 


submm position' 


counterpart position" 


infrared position 


Ar 


o 24 


5 100 


5160 


S 250 




S 500 




RA 


Dec 


RA 


Dec 


RA 


Dec 


n 


fljy 


mjy 


mJy 


mJy 


mJy 


mJy 


GN04 


12 36 16.60 


+62 15 20.00 S 


12 36 16.11 


+62 15 13.53 R 


12 36 16.10 


+62 15 13.58 


0.1 


302.8 ± 6.5 




12.5 ±1.9 


27.3 ±3.1 


25.7 ± 4.0 




GN05 


12 36 18.80 


+62 10 08.00 S 


12 36 19.13 


+62 10 04.32 R/M 


12 36 19.13 


+62 10 04.33 


0.0 


215.0 ±6.0 






9.3 ±3.1 






GN06 


12 36 18.70 


+62 15 53.00 S 


12 36 18.33 


+62 15 50.40 R 


12 36 18.33 


+62 15 50.41 


0.0 


330.0 + 7.6 


4.3 ± 1.0 


25.3 ± 2.0 


34.2 ±3.1 


46.7 ± 4.0 


27.4 ± 4.0 


GN07 


12 36 21.30 


+62 17 ll.OOS 


12 36 21.27 


+62 17 08.16 R 


12 36 20.98 


+62 17 09.54 


2.4 


366.9 ± 9.2 


4.3 ± 1.0 


10.5 ±2.0 


26.2 ±3.1 


27.6 ± 4.0 


13.2 ±3.9 


GN13 


12 36 50.50 


+62 13 17.00 S 


12 36 49.72 


+62 13 13.97 R/M 


12 36 49.72 


+62 13 12.88 


1.1 


371.0 ± 10.4 


9.0+1.9 


23.6 ± 2.4 


21.4±3.1 






GN15 


12 36 56.50 


+62 12 02.00 S 


12 36 55.82 


+62 12 01.13 M/1 


12 36 55.82 


+62 12 01.14 


0.0 


200.0 ± 6.0 






12.7 ±3.1 






GN19 


12 37 07.70 


+62 14 ll.OOS 


12 37 07.19 


+62 14 07.97 R 


12 37 07.19 


+62 14 07.98 


0.0 


276.7 ± 10.8 




10.1 ± 1.6 


21.8 ±3.1 


28.0 ± 4.0 


16.7 ± 3.9 


GN20 


12 3711.70 


+62 22 12.00 S 


12 37 11.88 


+62 22 12.11 R 


12 3711.88 


+62 22 12.09 


0.0 


68.9 ± 4.8 






12.6 ±3.1 






GN25 


12 36 28.70 


+62 10 47.00 S 


12 36 29.12 


+62 10 45.91 R 


12 36 29.11 


+62 10 45.91 


0.1 


724.0 ± 12.0 


7.6+1.0 


19.1 ± 1.9 


26.7 ± 3.0 


30.6 ±4.1 


12.4 ± 3.9 


GN26 


12 36 35.50 


+62 12 38.00 S 


12 36 34.51 


+62 12 40.93 R 


12 36 34.51 


+62 12 40.94 


0.0 


446.0 ±5.1 


34.4 ± 1.2 


66.0 ± 1.7 


54.3 ±3.1 


45.0 ± 4.2 


14.3 ± 4.3 


GN31 


12 36 53.10 


+62 1 1 20.00 S 


12 36 53.22 


+62 11 16.69 M 


12 36 53.22 


+62 11 16.70 


0.0 


367.0 ± 6.4 




5.7 ± 1.7 


13.9 + 3.1 


25.4 ± 4.0 


14.6 ± 4.0 


GN34 


12 37 06.50 


+62 21 12.00 S 


12 37 06.22 


+62 21 11.57 M 


12 37 06.22 


+62 21 11.60 


0.0 


81.7 ±4.1 




9.7 + 2.0 








GN20.2 


12 37 09.50 


+62 22 06.00 S 


12 37 08.77 


+62 22 01.78 R 


12 37 08.77 


+62 22 01.78 


0.0 


30.2 ± 5.6 




5.0 ± 1.0 


12.1 ±3.1 






GN39 


12 37 11.30 


+62 13 31.00 S 


12 37 11.33 


+62 13 31.02 R 


12 37 11.35 


+62 13 31.04 


0.1 


756.0 ± 9.3 


11.1 ±0.7 


33.5 ± 1.5 


53.5 ± 3.0 


39.9 ± 4.3 


25.9 ± 3.9 



Notes. 

' after the position, we indicate the (sub)mm observation from which the multi-wavelength identification has been made, S: SCUBA, L: Laboca, A: AzTEC and M: MAMBO. 
after the position, we indicate the nature of the observation that has provided the multi-wavelength identification of the (sub)mm source, R: radio, M: MIPS, I: IRAC. 



Table 4. submm, radio and redshift properties for our GOODS-N SMG sample. 













Submm properties 










Radio properties 


Redshift 


name 


SCUBA position 


5 850 


name 


AzTEC position 


S 1100 


name 


MAMBO position 


S 1200 


S 1.4GI1Z 






RA 


Dec 


mJy 




RA 


Dec 


mJy 




RA Dec 


mJy 






GN04 


12 36 16.60 


+62 15 20.00 


4.9 ± 0.7 


AzGN16 


12 36 16.18 


+62 15 18.10 


2.9+1.1 








89.5 ± 6.3 


2.578 


GN05 


12 36 18.80 


+62 10 08.00 


5.2+1.8 
















50.3 ± 15.4 


2.210 


GN06 


12 36 18.70 


+62 15 53.00 


7.5 + 0.9 


AzGN36 


12 36 17.38 


+62 15 45.50 


1.9+1.2 








178.9 ± 6.4 


1.865 


GN07 


12 3621.30 


+62 17 11.00 


8.9 ± 1.5 
















210.8 ± 12.9 


1.988 


GN13 


12 36 50.50 


+62 13 17.00 


1.9 ±0.4 
















45.4 ± 5.4 


0.475 


GN15 


12 36 56.50 


+62 12 02.00 


3.7 ± 0.4 
















16.2 ± 16.2 


2.743 


GN19 


12 37 07.70 


+6214 11.00 


8.0 ±3.1 
















77.0 ± 10.0 


2.484 


GN20 


12 37 11.70 


+62 22 12.00 


20.3 + 2.1 


AzGNOl 


12 37 12.04 


+62 22 11.50 


10.7 ± 0.9 


GN1200.1 


12 37 11.70 +62 22 11.00 


9.3 ± 0.5 


70.0+16.3 


4.055 


GN25 


12 36 28.70 


+62 10 47.00 


3.2 ± 1.4 
















93.8+ 12.9 


1.013 


GN26 


12 36 35.50 


+62 12 38.00 


2.2 ± 0.8 
















194.3 ± 10.4 


1.219 


GN31 


12 36 53.10 


+62 1 1 20.00 


2.1+0.6 
















16.3 + 5.4 


0.935 


GN34 


12 37 06.50 


+62 21 12.00 


3.8 ± 1.9 
















23.6 ± 23.6 


1.360 


GN20.2 


12 37 09.50 


+62 22 06.00 


9.9 ± 2.3 
















180.7 ± 8.4 


4.051 


GN39 


12 3711.30 


+62 13 31.00 


5.2 ± 2.4 


AzGN07 


12 37 11.94 


+62 13 30.10 


4.0 ± 1.1 


GN1200.3 


12 37 11.20 +62 13 28.00 


3.6 ± 0.6 


178.9 ± 8.6 


1.996 



Table 5. Mid- and far-infrared properties of our GOODS-S SMG sample. 



Reference submm source and its counterpart PEP/HerMES multi-wavelength counterpart 



name 


submm position'^ 


counterpart position" 


infrared position 


Ar 


524 




S 100 


S 160 


5250 


•S350 


5 500 




RA 


Dec 


RA 


Dec 


RA 


Dec 




yuJy 


mjy 


mJy 


mJy 


mJy 


mJy 


mJy 


LESSOlO 


03 32 19.02 


-27 52 19.40 L 


03 32 19.04 


-27 52 14.30 R 


03 32 19.05 


-27 52 14.46 


0.2 


161.9 ±4.6 






11.5 ±0.9 


24.0 ± 2.6 


30.3 ± 3.2 


20.8 ±4.3 


LESSOl 1 


03 32 13.58 


-27 56 02.50 L 


03 32 13.84 


-27 55 59.80 R 


03 32 13.85 


-27 55 59.93 


0.2 


119.1 ±6.0 




1.3 ±0.4 


3.4 ±0.6 


15.2 ±2.6 


19.2 ±3.4 


17.9 ±4.2 


LESS017 


03 32 07.59 


-27 51 23.00 L 


03 32 07.26 


-27 51 20.10R/M 


03 32 07.28 


-27 51 20.17 


0.3 


230.8 ± 5.6 




1.9 ±0.6 


7.4 ± 1.3 


25.0 ± 2.6 


29.0 ±3.2 


26.4 ± 4.2 


LESS018 


03 32 05.12 


-27 46 52.10 L 


03 32 04.87 


-27 46 47.40 R/M 


03 32 04.87 


-27 46 47.28 


0.1 


630.3 + 5.9 




3.7 ±0.5 


18.3 ±0.7 


40.6 + 2.6 


44.3 ± 3.6 


39.0 + 4.3 


LESS040 


03 32 46.74 


-27 51 20.90 L 


03 32 46.77 


-27 51 20.70 R/M 


03 32 46.78 


-27 51 20.90 


0.3 


I51.9±3.6 




2.2 ± 0.4 


7.1 ±0.9 


15.7 + 2.6 


1 9.0 ± 3.1 




LESS067 


03 32 43.28 


-27 55 17.90 L 


03 32 43.18 


-27 55 14.20 R/M 


03 32 43.19 


-27 55 14.36 


0.2 


583.5 ±5.1 


1.3 ±0.3 


4.4 ± 0.5 


15.0 ±0.8 


25.0 ± 2.6 


21.7 ±3.3 


24.9 ± 4.2 


LESS079 


03 32 21.25 


-27 56 23.50 L 


03 32 21.61 


-27 56 23.10 R/M 


03 32 21.60 


-27 56 23.34 


0.3 


615.6 ±7.0 




5.2 ± 0.6 


14.1 ±0.8 


33.9 ± 2.6 


30.1 ±3.4 


23.6 ± 4.3 



Notes. 

' after the position, we indicate the (sub)nim observation from which the multi-wavelength identification has been made, S: SCUBA, L: Laboca, A: AzTEC and M: MAMBO. 
after the position, we indicate the nature of the observation that has provided the multi-wavelength identification of the (sub)mm source, R: radio, M: MIPS. 



Table 6. submm, radio and redshift properties for our GOODS-S SMG sample. 









Submm properties 








Radio properties 


Redshift 


name 


LABOCA position 


5 870 


name 


AzTEC position 


5 1100 


5 1.4 GHz 






RA 


Dec 


mJy 




RA 


Dec 


mJy 






LESSOlO 


03 32 19.02 


-27 52 19.40 


9.1 ± 1.2 


AzTEC/GS2.1 


03 32 18.99 


-27 52 13.80 


6.3 ±0.5 


54.9 ±4.6 


2.437 


LESSOll 


03 32 13.58 


-27 56 02.50 


9.1 ±1.2 


AzTEC/GS7 


03 32 13.47 


-27 56 06.70 


3.5 ± 0.6 


55.1 ±4.2 


2.679 


LESS017 


03 32 07.59 


-27 5123.00 


7.6 ±1.3 


AzTEC/GSlO 


03 32 07.13 


-27 5125.30 


3.5 ± 0.6 


120.3 ± 4.3 


1.053 


LESS018 


03 32 05.12 


-27 46 52.10 


7.5 ±1.2 


AzTEC/GS8 


03 32 05.12 


-27 46 45.80 


3.1 ±0.5 


130.1 ±5.5 


2.214 


LESS040 


03 32 46.74 


-27 5120.90 


5.9 ± 1.3 


AzTEC/GS25 


03 32 46.96 


-27 5122.40 


1.8 ±0.5 


119.1 ±1.4 


1.593 


LESS067 


03 32 43.28 


-27 55 17.90 


5.2 ±1.4 










90.1 ±3.7 


2.122 


LESS079 


03 32 21.25 


-27 56 23.50 


4.7 ± 1.4 


AzTEC/GS23 


03 32 21.37 


-27 56 28.10 


2.1 ±0.6 


34.8 ±4.9 


2.073 



Table 7. Mid- and far-infrared properties of our LH SMG sample. 



Reference submm source and its counterpart PEP/HerMES multi-wavelength counterpart 



name 


submm position^ 


counterpart position" 


infrared position 


Ar 


o 24 


5 100 


S 160 


S250 




S350 


S 500 




RA 


Dec 


RA 


Dec 


RA 


Dec 


It 




mjy 


mJy 


mjy_ 




mJy 


mJy 


LOCK850.01 


10 52 01.42 


+57 24 43.04 S 


10 52 01.25 


+57 24 45.76 R 


10 52 01.25 


+57 24 45.90 


0.1 


178.5 ± 12.6 




15.2 + 3.1 


37.4 + 3 


9 


37.9 ± 5.4 




LOCK850.03 


10 52 38.25 


+57 24 36.54 S 


10 52 38.30 


+57 24 35.76 R 


10 52 38.30 


+57 24 35.77 


0.0 


175.0 + 23.0 


7.9 + 1.5 


19.3 + 2.5 










LOCK850.04 


10 52 04.17 


+57 26 58.85 S 


10 52 03.69 


+57 27 07.06 M 


10 52 03.70 


+57 27 07.40 


0.4 


1047.0 ± 17.0 


6.2 ± 1.0 


14.4 + 3.1 


40.0 + 1 


6 


50.7 ± 4.3 


29.6 ± 5.6 


LOCK850.12 


10 52 27.61 


+57 25 13.08 S 


10 52 27.58 


+57 25 12.46 R 


10 52 27.58 


+57 25 12.52 


0.1 


229.9 ± 10.5 


3.6 ± 1.0 


10.1 ±2.7 


22.7 + 3 


9 


28.2 + 5.5 




LOCK850.14 


10 52 30.11 


+57 22 15.55 S 


10 52 30.72 


+57 22 09.56 R/M 


10 52 30.74 


+57 22 09.62 


0.2 


194.1 ± 10.8 




10.4+ 1.8 


19.4 + 3 


8 


24.9 + 5.7 


30.6 + 8.3 


LOCK850.15 


10 53 19.20 


+57 21 10.64 S 


10 53 19.27 


+57 21 08.45 R 


10 53 19.28 


+57 21 08.67 


0.2 


380.8 ± 10.5 


5.9 ± 1.1 


11.4 + 3.1 


25.8 + 3 


8 


28.3 ± 6.7 




LOCK850.16 


10 51 51.45 


+57 26 37.00 S 


10 5151.69 


+57 26 36.09 R 


10 5151.68 


+57 26 36.08 


0.1 


336.6 ± 11.3 


15.1 + 1.3 


30.0 ± 1.7 


40.7 + 3 


8 


34.9 ± 5.4 




LOCK850.17 


10 51 58.25 


+57 18 00.81 S 


10 51 58.02 


+57 18 00.27 R 


10 51 58.03 


+57 18 00.29 


0.1 


256.3 ± 12.9 


9.1 ± 1.0 


28.4 ± 2.9 


35.7 + 3 


9 


30.5 ± 5.8 




LOCK850.33 


10 51 55.97 


+57 23 11.76 S 


10 5155.47 


+57 23 12.77 R 


10 51 55.46 


+57 23 12.89 


0.1 


112.9 ± 10.1 






18.6 + 3 


8 






SMMJ105238+571651 


10 52 38.19 


+57 16 51.10 S 






10 52 38.18 


+57 16 51.11 


0.1 


462.1 ± 10.6 


12.5 + 1.3 


15.0 + 2.9 


19.7 + 3 


8 






AzLOCK.l 


10 52 01.98 


+57 40 49.30 A 


10 52 01.92 


+57 40 51.50 R 


10 52 01.93 


+57 40 51.65 


0.2 


1473.5 ± 15.8 


12.2+1.5 


25.3 + 4.9 


59.8 + 3 


9 


56.4 + 6.2 


44.7 ± 9.4 


AzLOCK.5 


10 54 03.76 


+57 25 53.70 A 


10 54 03.75 


+57 25 53.50 R 


10 54 03.79 


+57 25 53.62 


0.3 


225.9 ± 20.3 


3.9+1.1 


31.2 + 2.7 


57.3 + 3 


9 


51.1 ±5.9 




AzLOCK.lO 


10 54 06.44 


+57 33 09.60 A 


10 54 06.83 


+57 33 09. 10 R 


10 54 06.86 


+57 33 09.42 


0.4 


622.8 ± 16.4 


4.5 + 1.3 


9.4 + 2.8 


12.4 + 3 


9 


34.5 ± 7.5 


35.0 ± 8.7 


AzLOCK.62 


10 52 11.61 


+57 35 10.70 A 


10 52 11.85 


+57 35 10.50 M 


10 52 11.85 


+57 35 10.49 


0.0 


243.2 ± 11.0 




14.2 ± 2.6 


18.1 + 3 


8 


26.5 ± 5.6 





Notes. 

' after the position, we indicate the (sub)nun observation from which the multi-wavelength identification has been made, S: SCUBA, L: Laboca, A: AzTEC and M: MAMBO. 
after the position, we indicate the nature of the observation that has provided the multi-wavelength identification of the (sub)mm source, R: radio, M: MIPS. 



Table 8. submm, radio and redshift properties for our LH SMG sample. 













Submm properties 














Radio properties 


Redshift 


name 


SCUBA position 




name 


AzTEC position 


1100 


name 


MAMBO position 


S 1200 


S 1.4 GHz 






RA 


Dec 


mJy 




RA 


Dec 


mJy 




RA 


Dec 


mJy 






LOCK850.01 


10 52 01.42 


+57 24 43.04 


8.8+1.0 


AzLOCK.8 


10 52 01.14 


+57 24 43.00 


4.7 + 1.0 


LE1200.5 


10 52 01.30 


+57 24 48.00 


3.4 + 0.6 


110.0 + 6.0 


3.380 


LOCK850.03 


10 52 38.25 


+57 24 36.54 


10.9+ 1.8 


AzLOCK.24 


10 52 38.46 


+57 24 36.80 


3.0+1.0 


LEI 200.1 


10 52 38.30 


+57 24 37.00 


4.8 + 0.6 


25.8 + 4.9 


3.036 


LOCK850.04 


10 52 04.17 


+57 26 58.85 


10.6+1.7 


AzLOCK.7 


10 52 03.89 


+57 27 00.50 


4.8 + 0.9 


LE1200.3 


10 52 04.10 


+57 26 58.00 


3.6 + 0.6 


47.0 + 5.7 


1.480 


LOCK850.12 


10 52 27.61 


+57 25 13.08 


6.1 + 1.7 










LE1200.6 


10 52 27.50 


+57 25 15.00 


2.8 + 0.5 


44.3 + 5.1 


2.470 


LOCK850.14 


10 52 30.11 


+57 22 15.55 


7.2 + 1.8 










LE1200.10 


10 52 29.90 


+57 22 05.00 


2.9 + 0.7 


37.4 + 4.2 


2.611 


LOCK850.15 


10 53 19.20 


+57 21 10.64 


13.2 + 4.3 


AzLOCK.17 


10 53 19.47 


+57 21 05.30 


3.6+1.0 










105.4 + 5.0 


2.760 


LOCK850.16 


10 5151.45 


+57 26 37.00 


5.8 + 1.8 


















106.0 + 6.0 


1.620 


LOCK850.17 


10 51 58.25 


+57 18 00.81 


4.7 + 1.3 










LE1200.11 


10 51 58.30 


+57 17 53.00 


2.9 + 0.7 


92.3 + 4.5 


2.694 


LOCK850.33 


10 51 55.97 


+57 23 11.76 


3.8 + 1.0 










LE1200.12 


10 5155.50 


+57 23 10.00 


3.3 + 0.8 


51.0 + 4.3 


2.686 


SMMJ105238+571651 


10 52 38.19 


+57 16 51.10 


5.3 + 1.6 


















71.1 + 12.6 


1.852 










AzLOCK. 1 


10 52 01.98 


+57 40 49.30 


6.6 + 0.9 










258.0 ± 11.0 


2.500 










AzLOCK.5 


10 54 03.76 


+57 25 53.70 


4.9+1.0 










138.0 + 9.0 


2.820 










AzLOCK. 10 


10 54 06.44 


+57 33 09.60 


4.1 + 0.9 










77.0 + 9.0 


2.560 










AzLOCK.62 


10 52 11.61 


+57 35 10.70 


2.0+1.0 












2.480 



Table 9. Mid- and far-infrared properties of our COSMOS SMG sample. 



Reference submm source and its counterpart PEP/HerMES multi-wavelength counterpait 



name 


submm position" 


counteipart position^ 


infrared position 


A/' 


■> 24 


^ 100 


^ 160 


■J 250 


^ 350 


^ 500 




RA 


Dec 


RA 


Dec 


RA 


Dec 


If 




mJy 


mJy 


mJy 


mJy 


mJy 


COSLA-121R1I 


09 59 38.82 


+02 08 41.28 L 


09 59 38.94 


+02 08 49.49 R 


09 59 38.96 


+02 08 48.98 


0.6 


325.0+ 16.0 






18.0 + 2.7 


17.5 ±3.3 




COSLA-127R1I 


10 01 24.58 


+01 56 06.93 L 


10 01 23.86 


+01 56 13.39 R 


10 01 23.88 


+01 56 13.49 


0.4 




11. 2± 1.4 


33.1 +4.4 


28.8 + 2.7 


33.4 + 4.1 




COSLA-155R1K 


09 59 39.06 


+02 21 21.19L 


09 59 39.06 


+02 21 26.42 R 


09 59 39.06 


+02 21 26.53 


0.1 








11.2 + 2.7 


22.7 + 3.6 




COSLA-I63RII 


09 59 28.12 


+02 07 48.45 L 


09 59 28.63 


+02 07 49.49 R 


09 59 28.57 


+02 07 48.07 


1.7 


184.0 + 35.0 


6.9+1.6 




22.4 + 2.7 


23.6 + 5.1 


33.5 + 5.0 


COSLA-012RII 


10 00 30.16 


+02 41 37.61 L 


10 00 30.25 


+02 41 46.35 R 


10 00 30.33 


+02 41 46.46 


1.2 


1095.0+ 163.0 


17.6+1.8 


40.4 + 3.4 


65.2 + 2.7 


72.0 + 3.8 


39.8 + 6.0 


AzTECJ100008+022612 


10 00 07.95 


+02 26 08.16 A 


10 00 08.05 


+02 26 12.20 S 


10 00 08.10 


+02 26 11.57 


1.0 


287.0+ 15.0 




11.8 + 3.5 


21.3 + 2.7 


37.0 + 4.5 




AzTECJI00019+023206 


10 00 19.75 


+02 32 04.40 A 


10 00 19.77 


+02 32 04.33 S 


10 00 19.77 


+02 32 03.95 


0.4 


189.0+13.0 




31.4 + 3.7 


31.7 + 2.7 


32.9 ±4.8 




AzTECJ100020+023518 


10 00 20.70 


+02 35 20.50 A 


10 00 20.70 


+02 35 20.50 S 


10 00 20.95 


+02 35 18.80 


3.9 








16.9 + 2.7 


17.8 ±3.5 




AzTECJ 1 00008+024008 


10 00 08.91 


+02 40 09.60 A 


10 00 08.94 


+02 40 10.70 S 


10 00 08.95 


+02 40 10.67 


0.2 


660.0+ 17.0 


10.5 + 1.7 


45.1+3.9 


77.9 + 2.7 


73.3 + 3.3 


52.0 + 5.0 


AzTECJ095939+023408 


09 59 39.30 


+02 34 08.00 A 


09 59 39.18 


+02 34 03.67 R 


09 59 39.20 


+02 34 02.75 


0.9 








29.2 + 2.7 


28.3 + 3.3 


20.9 + 5.2 


MAMBOll 


10 00 38.10 


+2 08 25.00 M 


10 00 38.01 


+02 08 22.57 R 


10 00 38.01 


+02 08 22.49 


0.1 


1589.0+ 115.0 


27.4+ 1.7 


72.9 + 5.1 


81.9 + 2.7 


69.3 + 3.3 


38.7 + 4.8 



Notes. 

" after the position, we indicate the (sub)mm observation from which the multi-wavelength identification has been made, S: SCUBA, L: Laboca, A: AzTEC and M: MAMBO. 
^ after the position, we indicate the nature of the observation that has provided the multi-wavelength identification of the (sub)mm source, R: radio, S: SMA. 



Table 10. submm, radio and redshift properties for our COSMOS SMG sample. 



index 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



Submm properties 



name 


LABOCA position 


^ 870//m 




RA 


Dec 


mJy 


C0SLA121 


09 59 38.82 


+02 08 41.28 


11.8 + 3.8 


COSLA127 


10 01 24.58 


+01 56 06.93 


15.8 + 5.1 


COSLA155 


09 59 39.06 


+0221 21.19 


13.0 + 4.3 


COSLA163 


09 59 28.12 


+02 07 48.45 


13.3 + 4.4 


C0SLA12 


10 00 30.16 


+02 41 37.61 


23.4 + 5.6 


C0SLA2 


10 00 07.95 


+02 26 08.16 


18.6 + 3.9 



SMA position 
RA Dec 



■^890 

mJy 



name 



RA 



AzTEC position 
Dec 



SMA2 
SMA5 
SMA3 
SMAll 



10 00 08.05 
10 00 19.77 
10 00 20.70 
10 00 08.94 



+02 26 12.20 
+02 32 04.33 
+02 35 20.50 
+02 40 10.70 



12.4+1.0 
9.3 + 1.3 
8.7+1.5 
14.4 + 4.2 



AzTECJ 100008+0226 12 
AzTECJ 1000 19+023206 
AzTECJ100020+023518 
AzTECJ 1 00008+024008 
AzTECJ095939+023408 



10 00 08.05 
10 00 19.75 
10 00 20.70 
10 00 08.91 
09 59 39.30 



+02 26 12.20 
+02 32 04.40 
+02 35 20.50 
+02 40 09.60 
+02 34 08.00 



^ 1100 

mJy 



AzTECJ095939+022124 09 59 39.01 +02 21 24.50 1.3 + 0.5 



8.3 + 1.3 
6.5 + 1.2 
5.9+1.3 

4.7 + 1.3 

3.8 + 1.4 



I. 



o 



Table fTOl continued. 



Submm propeities Radio propeities Redshift 

index name MAMBO position i^i200 '^i.4GHz 

RA Dec mJy yujy 



1 ... ... 54.0+ 10.0 1.850 

2 ... ... 51.0+ 10.0 0.907 

3 ... ... 128.0 + 51.0 2.974 

4 ... ... 54.0+10.0 1.178 

5 ... ... ... 1.260 

6 ... ... ... 1.120 

7 ... ... 140.0 + 30.0 3.971 

8 ... ... ... 5.310 

9 ... ... 140.0+ 12.0 1.599 

10 ... ... 69.0+ 11.0 0.834 

11 MAMBOll 1000 38.10 +02 08 25.00 4.6 + 0.9 237.0 + 27.0 1.830 



Table 11. Mid- and far-infrared properties of our lensed-SMG sample. 





Reference subnun source and its counterpart 












PEP/HerMES multi-wavelength counterpart 








field 


name 


submm position" 


counterpart position 


infrared position 


Ar 


J 24 


5 100 


5 160 


•5250 


•S350 


•Ssoo 






RA 


Dec 


RA 


Dec 


RA 


Dec 


It 




mjy 


mJy 


mJy 


mJy 


mJy 


A1835 


SMMJ14011+0252 


14 01 04.96 


+02 52 23.50 S 


14 01 04.97 


+02 52 24.60 


14 01 04.97 


+02 52 24.54 


0.1 


883.4 ± 14.5 


11.6 ±0.8 


33.5 ± 1.4 


61.7 + 6.0 


63.1 ±3.8 


48.5 ± 4.0 


A1835 


SMMJ 14009+0252 


14 00 57.55 


+02 52 48.60 S 


14 00 57.57 


+02 52 49.10 


14 00 57.58 


+02 52 48.95 


0.2 


297.7 ± 13.7 


4.9 ± 0.7 


27.7+1.9 


66.4 ± 6.0 


65.8 + 3.8 


53.7 ± 4.2 


A2219 


SMM J 16403+4644 


16 40 19.40 


+46 44 Ol.OOS 


16 40 19.50 


+46 44 00.50 


16 40 19.46 


+46 44 00.81 


0.5 


729.7 ± 7.6 


17.1 ±0.9 


35.6+1.9 


50.5 + 5.8 


44.9 + 3.6 


34.7 ± 6.5 


MS 1054 


SMMJ 10570-0334 


10 57 02.20 


-03 36 04.00 S 






10 57 02.50 


-03 36 02.52 


4.7 


253.1 + 12.0 




10.6 + 2.1 


15.2+1.2 


17.0+ 1.1 




CL0024 


SMMJ00266+1708 


00 26 34.10 


+ 17 08 32.00 S 


00 26 34.06 


+ 17 08 33.10 


00 26 34.09 


+ 17 08 34.16 


1.1 


283.4+ 15.4 


4.7 ± 0.7 


24.5 ± 2.2 


52.8 + 7.2 


61.3 + 4.5 


44.6 + 4.7 


MS0451 


SMMJ04542-0301 


04 54 12.50 


-03 01 04.00 S 






04 54 12.73 


-03 01 09.25 




169.9 ± 10.5 


14.6 ± 1.2 


40.6 ± 7.7 


76.0 ± 15.1 


94.4 ± 8.8 


84.1 ± 11.5 


A2390 


SMMJ21536+1742 


21 53 38.20 


-H17 42 16.00 S 


21 53 38.35 


+17 42 20.70 


21 53 38.51 


+17 42 17.73 


3.8 


382.7 ± 8.0 


3.8 ± 0.9 


11.9 + 2.1 


52.4 ± 12.7 


54.8 ±7.1 


47.5 ± 7.2 


A2218 


SMMJ16354+6611 


16 35 41,20 


+66 11 44.00 S 


16 35 41.20 


+661144.00 


16 35 40.73 


+661142.96 


3.1 


411.8 ± 156.3 


1.0 ± 0.7 


13.6 + 1.2 


20.5 ± 6.0 


25.7 ±3.0 


17.4 ±3.3 


A2218 


SMMJ16355+66120C 


16 35 50.96 


+66 12 05.50 S 


16 35 50.96 


+66 12 05.50 


16 35 50.75 


+66 12 06.31 


1.5 


280.2 ± 26.2 


4.9 ± 0.7 


18.5 ± 1.6 


37.8 ± 6.0 


35.9 ±3.1 




A2218 


SMMJ16355+66122B 


16 35 54.10 


+66 12 23.80 S 


16 35 54.10 


+66 12 23.80 


16 35 53.91 


+66 12 22.43 


1.8 


2709.2 ± 129.0 


9.4 + 0.9 


28.1 + 1.9 


78.1 + 6.0 


70.2 ± 7.2 


74.7 ± 3.5 


A2218 


SMMJ16355+66123A 


16 35 55.18 


+66 12 37.20 S 


16 35 55.18 


+66 12 37.20 


16 35 54.96 


+66 12 37.04 


1.4 


2955.0+ 123.4 


5.9+1.1 


21.3 + 2.0 


53.7 + 6.0 


47.1 + 3.1 


26.9 + 3.5 


A2218 


SMMJ16355+6611 


16 35 55.20 


+66 11 50.00 S 


16 35 55.20 


+66 1 1 50.00 


16 35 55.00 


+66 1 1 50.62 


1.4 


4724.5 ± 182.7 


26.4 + 0.8 


51.0+1.5 


42.1 +6.0 


37.2 + 3.0 


14.9 + 3.3 


A370 


SMMJ02399-0136 


02 39 51.90 


-01 35 59.00 S 


02 39 51.88 


-01 35 58.00 


02 39 51.86 


-01 35 58.45 


0.5 


1231.0+ 11.6 


13.0 + 0.8 


30.1 ± 1.8 


65.7 + 6.1 


72.0 ± 3.9 


62.8 + 4.6 


A370 


SMMJ02399-0134 


02 39 56.40 


-01 34 27.00 S 


02 39 56.51 


-01 34 27.10 


02 39 56.58 


-01 34 26.10 


1.4 


2402.0 ± 10.9 


45.8 ± 0.8 


111.1 ± 1.9 


125.8 + 6.1 


89.8 + 3.9 


41.7 + 4.4 


A1689 


SMMJ131 15-1208 


13 11 29.10 


-01 20 49.00 S 


13 11 29.14 


-01 20 46.50 


13 11 29.14 


-01 20 46.47 


0.1 


329.8 ± 12.8 


3.3 ± 0.4 


7.6 + 0.8 




16.1 ±4.1 


17.3 + 4.4 



Notes. 

^ after the position, we indicate the (sub)mni observation from which the multi-wavelength identification has been made, S: SCUBA, L: Laboca, A: AzTEC and M: MAMBO. 



Table 12. submm, radio, redshift and magnification properties for our lensed-SMG sample. 







Submm properties 








Radio properties 


Redshift 


Magnification 


field 


name 


SCUBA position 




■^850 


^ 1350 


•5 1.4GHz 










RA 


Dec 


mJy 


mJy 


mJy 


/zJy 






A1835 


SMMJ14011+0252 


14 01 04.96 


+02 52 23.50 


41.9 + 6.9 


14.6 ±1.8 


6.1 ± 1.5 


115.0 ±30.0 


2.565 


3.00 


A1835 


SMMJ 14009+0252 


14 00 57.55 


+02 52 48.60 


32.7 ± 8.9 


15.6+1.9 


5.6 ± 1.7 


529.0 ± 30.0 


2.934 


1.50 


A2219 


SMMJ 16403+4644 


16 40 19.40 


+46 44 01.00 


53.4+16.0 


10.0 + 2.0 


4.3 ± 2.0 




2.030 


3.60 


MS1054 


SMMJ 10570-0334 


10 57 02.20 


-03 36 04.00 


10.9+ 10.9 


4.4+1.4 






2.423 


1.10 


CL0024 


SMMJ00266+1708 


00 26 34.10 


+ 17 08 32.00 




18.6+1.5 




94.0 ± 15.0 


2.730 


2.40 


MS0451 


SMMJ04542-0301 


04 54 12.50 


-03 01 04.00 




13.8 ±2.0 






2.911 


50.00 


A2390 


SMMJ21536+1742 


21 53 38.20 


+17 42 16.00 




6.7 ± 2.2 






1.020 


1.90 


A2218 


SMMJ16354+66114 


16 35 41.20 


+66 11 44.00 


53.4 ± 16.0 


10.4 ± 1.4 






3.188 


1.70 


A2218 


SMMJ16355+66120C 


16 35 50.96 


+66 12 05.50 


22.9 ± 6.9 


8.7 ± 1.1 






2.516 


9.00 


A2218 


SMMJ16355+66122B 


16 35 54.10 


+66 12 23.80 


46.4 ± 13.9 


16.1 ± 1.6 






2.516 


22.00 


A2218 


SMMJ16355+66123A 


16 35 55.18 


+66 12 37.20 


31.8 + 9.5 


12.8+1.5 






2.516 


14.00 


A2218 


SMMJ16355+6611 


16 35 55.20 


+66 1 1 50.00 


17.1 ± 5.1 


3.1+0.7 






1.034 


7.60 


A370 


SMMJ02399-0136 


02 39 51.90 


-01 35 59.00 


85.0+ 10.0 


23.0 + 2.0 




526.0 ± 10.0 


2.810 


2.50 


A370 


SMMJ02399-0134 


02 39 56.40 


-01 34 27.00 


42.0+ 10.0 


11.0 + 2.0 




500.0 ± 10.0 


1.060 


2.50 


A1689 


SMMJ131 15-1208 


13 11 29.10 


-01 20 49.00 


21.0 + 6.4 


4.7 + 0.8 






2.630 


21.60 
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Table 13. Dust properties of our SMGs. 



Multi-r'' 



Field 


Name 




logCLiR)"" 


Tc 


log(Mdusl) 


log (Lir)' 






K 


Lo 


K 


Mo 


Lo 


GOODSN 


GN04 


41 ± 1 


12.91 ±0.11 


27+ 1 


8.40 ± 0.05 


12.80 ±0.07 


GOODSN 


GN05 


27 ±3 


12.20 + 0.49 


19+1 


8.75 ± 0.25 


12.23 ±0.30 


GOODSN 


GN06 


35 ± 1 


12.77 + 0.06 


23 + 1 


8.75 ± 0.05 


12.65 ±0.05 
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Notes. 

" In this model the dust emissivity index is fixed to /? = 1.5 (see text for details). 
^ In this model we use yff - 2.0, y - 7.3 and R = 3 kpc (see text for details). 

The infrared luminosities of our lensed-SMGs have been de-magnified using the magnification factors given in Table |12| 
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Fig. A.l. Spectral energy distribution of our SMGs. Red diamonds present tlie PACS and SPIRE measurements, wiiile green squares present 
multi-wavelengtli ancillary data taken from the literature. The modified blackbody emission (fl = 1.5) best-fitting the data are shown by dashed 
blue lines. The power-law temperature distribution model (J3 = 2.0, y = 7.3 and R = 3kpc) which best-fits the data are shown by solid red lines. 
Dotted lines present the CEOl SED template which best-fits the far-infrared observations. 
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. A.4. Figure lATTI continued. 



Magnelli et al.: Far-infrared properties of SMGs 



10" 



10" 



10'' 

10^ 

10' 
lO'' 

10' 



lo-* 

10^ 

10' 
10^ 

10' 



SMMJ16403+46440 



z = 2.0 
T =25 K 



L,R=2.6 xio''' Lq 

X^multi-T = 6.7 




SMMJ105702. 2-033404 



z = 2.4 

T =25 K 



X inulli-T^2.3 




SMMJ00266+1708 
z = 2.7 
T =25 K 



X multi-T~^-^ 




SMMJ04542-0301 

z = 2.9 LiR=5.7 xio" Lq 

T,=25 K X'mum-T=34.1 



SMMJ21536+1742 
z=1.0 
T =14 K 



SMMJ163541. 2 + 661144 





--38.0 



z = 3.2 
T =27 K 



= 5.8 xlO''^ Lq 

X'multi-T=20.5 




SMMJ 163550. 9 + 661207C 

z = 2.5 L,R=9.5 xio" La 



T =23 K 



X multi-T 42.7 



SMMJ163554.2 + 661225B 



z = 2.5 
T =21 K 



L,R=7.6 XlO" L;; 

X inulti-T='''3- 



SMMJ163555.2+661238A 



z = 2.5 
T =22 K 



L,R=8.1 xio" Lq 

X multi-T=25. 1 



lo-' 



10' 

10' 
10'^ 

10' 




4- 

SMMJ163555. 2 + 661150 





z=1.0 
T =21 K 



L,„=3.5 xlO" L„ 



L„=8.7 xlO' 



L,o=2.0 xlO'^ Ln 



X multi-T" 



= 16.7 



X mulU-T=36.3 



X multi-T~ 17.2 



10'' 
10^ 
10' 




10' 



lO"" 



10'' 



10^ 



10' 



10' 



[/j.m] 



lO'' 



10* 



10' 



10' 10"^ 10^ 10* 10' 

observed 



. A.5. Figure lATTI continued. 
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